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FOREWORD 


$ 


This thermal network modelling handbook was generated-under Task Ho. 4 
efnNASA Contract HAS 9-10435 entitled "Advanced SINQA Thermal Analyz 

Development." / 

Under the.-total coniract^three tasks were devoted w "" 
Task Ho. 1 for: development of "the User's_Manual, TRW Report"Ho.l46SQ.-H001 
Xoo. April 1971 * -Task No. 2 -for development.of the Program Manual . TO 
Report Ho- 14690-H002-RD-00, June 1321. and. Task No- 4_for dove 
•i.e Thermal Network Modelling Handbook contained herein. 

Contrasted to the User's Manual and Program Manual which describe the 
complete usage and content of the SINDA Computer Program, a c ^ 1 * t V 
manuscript describing. all the comities of Thermar M^he^ lcal^dell 

was far beyond the scope and budgetary alloeatlons of this ta.k. Hence, 
thi s^ document was 1nte“ded-to serve as an Initial building block for a 
future comprehensive document. Using this guideline a three-fo d purpose 
r«t ZLd: Cl) to acquaint the new user with the terminology and 

c co e in thermal mathematical modeling. (2) to P^ent th^re 
experienced and occasional user with ,uUk formulae and metho^s or solving 
everyd«tproblems, CeupHd^lth study cases which lend Insfght Into the 
relationships, that- exist ,™n,_th^varlous solution technlquesandparam- 
eters, andU) to begin- to. catalogue In an orderly, fashion^ comnmn 
formulae whtch-ln- the futu re m be ap plied to automated conversational 

1 anguage techniques. 




m 


TABLE OF CONTENTS 


J 


Page 

LO INTRODUCTION T 

1.1 Thermal Math -Model 1 ng as a Cognitive Process . . . . « l~ 

2.0 THERMAL MATHEMATICAL MODELING . . 3 

2.1 Network Solution-—. .... ..... . . . . 3 

2.2- Modeling Elements^. . 4- 

2.2-. 1 ...Jtodes..-..*-*-* . . 4 . ...4 . • — . . 4 - 

2. 2. 1.1 Concepts . 4 

2.2.1.2 Node Types . ..... 5 

2.2-.1.3 Method of Nodal! nation 7 

2. 2. 1.4 Computational Methods-Nodes 11 

2.2.2 Conductors. . . . . — . — 12 

2. 2. 2.1 Concepts ... 12 

2. 2JJ. 2 Conductor Types — 12 

2. 2. 2. 3 Computational Methods - Conduction 

Conductors 15 

2.2.2.3J Rectangular Nodes . — . . . *5 

2.2. 2- 3*2— Circular Sections . — , 16 

2. 2.2. 3-3- Parallel Conductors . . .. - 16. 

2. 2.2- 214— Series- Conductors . « . . 17 

2. 2. 2. 4 Computational Methods - Convection- 

Conductors ..... ...... ; 18" 

2t 2, 2.4.1 Combined Natural^and ... 

ForeecL Convection— .... 19 

2 ~. 2. 2. 4. 2 Natural- Convection 

Equations . , , r-. . r . 2+ - 

2. 2.2. 4. 3 Forced Convection in Tubes- - 

and Suctsr-r .... j.- . . 26- 

2. 2. 2. 4. -4 Forced Convec-tlon-Over 

Flat Pl ates - 31 


) 



1 v 


i 


/ 




Vj' 




Page 

2.2. 2.4.5— Fo* :.5d Convection 

r . Cylinders. 33_ 

2.2. 2. 4; 6 rorced Convection 

Over Spheres-.^ 36 

2.2.2JS" Computational-Methods - Radiation 

Condaeto rs. .......... ... . . „ . „ 38 

2.2. 2. 6 Computational Methods —Mass. flow 

Conductors—. . i , . . , ^ 57 


2.2.3. Energy Sources or. Sinks — . . . 5ft 

2. 2. 3.1 Concepts 58 

2. 2. 3. 2 Types of Heat Sources or Sinks 58 

2. 2. 3. 3 Computational Considerations - 

Sources or Sinks 58— 

2.3 Network Solution— . . 59 

2.3.1 Steady State . 61 

2.3.2 Transient Analysis . 63 

2.3.2.1 - Forward Differencing. . , 63 

2. 3.2 12. Backward Differencing f>5 

2.3.3 Sunmary^ of "Other Techniques. . . . 67 

2JL. Model ing-Parameters- fiR 

3.0 OPERATIONAL PARAMETER RELATIONSHIPS . 70. 

3.1 One-JBimensional Bar of Metal, w 70 

3.2 Other 0ne=d)1aiens1onal Cases-. . . . — . . . * . .- . . 78 

3.3 Two-Dimensional Plate of=Meta1. . * . — . . * 78 

3.4 Other Two-Dimensional Cases ; . . . . 79 

4.0 REFERENCES . . . . 82 



v 


LIST OF FIGURES 


2=7 
2-8 
2-9 — 
2-10 
2-U 
2-12 
2-13 
2-14 
2-15 
2-16 
2-17 


2-19 

2-20 

2-21 

2-22 

2^23 

2=24 


Figure 

2-1 Nodalization. . ....»•*• . 

2=2 Temperature- Distributions 

2=3 Alternate- Nasalization Methods. 

2=4 Polygonal. Nodal 1 zation vs Rectangular Nodalization 
2-5 Nodalization of Circular Elements-. ... 

2-6. Sample. Boundary Node* ... * »-*•*• * 

Use of Arithmetic. Nodes t o Mode l. Surfaces 
Thermal Network Elements. .... * ...... ♦ * • - 

Conduction Conductor. .......... 

Convection Conductor 

Radiation Conductor . 

Mass Flow Conductors. » • ~ • 

Simple Conductor Representing Heat Flow Path 

Area and Length-Equivalents for Circular Section Nodes. 

parallel Condirctor Flow Paths 

Series Conductor Flow-Paths ....... 

Empirical Constants for Natural Convection - 

Laminar-Flow. •— -* 

2-4 a— Empirical Constants for Natural Convection - 
Turbulent Flew. . . ... - * . . . * » • * - 
Empirical' Constants for Forced Convection Transition 
Flow in Tubes and- Ducts' *— *- » » *- * • *— 

Flow -Over Cylinders . . . . * . . • • 

Empirical Constant-for Cylinder Stagnation-Point. . . 

Radiant Interchange-Configuration Factors — 

Point Sources 

Radiant Interchange Configuration Faetors - 

Line-Sources. ♦ ; . ; 

Radiant- Interchange- Configuration Factors - - 

Plane Sources * 

2-25 Forward Extrapolation- Used in Forward Differencing.. . 
2-26 Backward- Extrapolation Used in- Backward Differencing. 


Page 

4 

5 
8 


-a .* 


29 

33- 

34 

40 

46 

48 

63 

65 


„_8‘ 


. . 

9 



9 

10 



12 


•*. 

12 

13 



14 



14 


J 

15 


i 

16 

h 

r> 

16 

17 


24 


i 

■ -■£ 

25" 




Vi 



Figure 

3-1 

3-2 


3-3 

3*4 

3-5 


3-6- 

3-7 


3-8 


Page 


Bar of Mf*taT. - • •••♦•*•• 

The Effect of a. Variation- In Size of Time Step for 
TwsTveJtode Model at t?25 Seconds and- x*l Inch. _ . . . . 

The Effecto-Pe, Variation. In Node Size at-t~25 Seconds 
and Inch.. . . » * < * * • — - *-«-—* — • — * * • r * * * * 

The Effect, of-a Variation In- Size of^Ttme-Step for 
Twelve_Node_Model atit»100 Seconds and x=6 Inches . . . — » 

The^Effect- of' a. Variation Jn Node_Stze_at t®T0Q.'_ 

Seconds and' x»6 Inches . . . < . . * • ■ *-• * 

Boundary Conditions on. Two-Dimensional Plate* 

The Effect of Variation in Time-Step for Center of Plate 

at t-tflO Seconds^ 

The Effect of a Variation, in Node Size for- Center ef 
Plate at t«10Q Seconds * 


.70 

74 

75 

76 

71 

7ft 

80 

81 




vft 


LIST OF TABLES 


Table Pa 9 e 

2-Y_ Thermal=Electri cal System Analogy. 3 

Z-2 OrderofMagnitude-cf Convective Heat Transfer- 

Coefficient* . . .■ . •— *. * — • — *• • • » • i® 

2-3 Empirical Constants* for Lamina r BowJJver aCylimler--* . 35 

2-4 Effective Emlttance for Flat Platestlnside-a Black 

EnclDSUFe*...*..*..*,..* 51 


viti 


NOMENCLATURE 


A 

a 

Area 

C— 

_B 

Thermal Capacitance 

Cp 

3 

Specific-heat 

D 

a 

D1 ameter 

dh 

m 

Hydraulic diameter 

r 

m. 

Voltage 

F 

= 

Radiation, configuration, (form) factor 

G 

* 

ThermaT~conductance . 

h 

9 

Convective heat transfer coefficient- 

L 

m 

Length or running length 

I 

B 

Current 

k 

9 

Thermal conductivity 

rfi 

9 

Mass flow rate 

n 

s 

Arbitrary- exponent- 

N 

3 

Number of Iterations 

* 

a 

‘ Meat- rate 

r 

s 

Radi us 

R — 

c 

Resistance 

t 

m 

Time- 

T 

a 

Temperature 

Too - 

s 

-Surrounding media or free- stream temperature- 

U- 

-a* 

Vetocfty 

Uoo 

a 

Free stream velocity 

V- 

«*- 

Vo lime 

# 

» — FI ow -rate 

Ws_ 

a 

Sampling frequency- 

W C 

a 

Maximum frequency component 

X 

a 

Arbitrary distance- 
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NOMENCLATURE (Concluded) 


aT Temperature difference 
§ = Script-F (grey beefy from factor) 

o *_ -Thermal, dlffusivlty 
6 » Coefficient of' Volumetric- expansion 

p * Density 

a * Stephan-Bol tzmann constant* 
t * Stability- factor 

e = Emittance 

©- ® Angle 

x * Radi ati on. 1 1 neari zatl on factor 

6 * Convergence criterion-^ relaxation criterion) 

c = Damping factor 

Symbols, subscripts and units not specifically mentioned In the 
Nomenclature are explained at the point of usage-within the text. 
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1.0 INTRODUCTION 


1 .1 Thermal Hath Model 1ng~as a Cognitive Process 

A brief introduction to the_rudimentaify techniques of thermal modeling 
coup led .with a -simple understanding of the various, basin heat transfer 
mechanisms are-the prerequisites fora beginning -thermal math, modeler. 
Properly applied*, the -body of concepts, principles and techniques~applicable 
to. thermal math. modeling constitutes a valid enqineeHng- tool. which can be 
applied. to the solution of real engineering problems. A good lumped param- 
eter representation of a_thermal system requireswin addition to the basic 
principles' and techniques, an- elusive mixture of experience (with real 
systems^, both physical and model) and Engineering. Judgement to transfer the 
end product into an accurate, versatile and cost effective thermal math., 
model . 

Generally, the problems encountered in developing a thermal. math model- 
reduce-Jto an overall object Of achieving the greatest accuracy for the- least 
cost. Cost factors are rather well defined, and fall into two classes: 

(1) the cost of developing the model, and (2) the cost of using the model 
Development costs ean be based almost solely on the actual engineering man- 
power required to do the_job within, the constraints of time- and budget. 
However, the potential costs Involved in. using a model are often not as 
obvious-nor as- Uneat— For^xample, most thermal math-models will bt ,ne 
ly?ed on a computer which is highly prioritized^. Ill such an environment, 
a-celatl vely- small increase In required processor time,-perhaps-from Sk. 
mlnutes— to 15 minutes, often results In a reduced job priority and a corres- 
ponding slowdown in turnaround time, perhap&_from One dqy to one- week. 

The- problem of - achieving accuracy in a math model , while subjeefc to 
coot constraints, varies greatly from- one- thermal math model to another. 
General accuracy requirements may be-as straight forward as "temperature 
accuracy shalV be-compattble-With- thermocouple A/D converter quantization 
error," or thermostat hysteresis." 0n the Other hand#- accuracy levels 


might be Indirectly indicated by^requirtng that a model must "be sufficient- 
ly detailed to permit meaningful paraaietrlc analyses-wtth respect to Insu- 
lation thickness variations in Increments of 1/4 Inch," Clearly, thet£..ts 
going to be a great- deal of Engineering. Judgement involved In .developing a 
model that is “sufficiently detailed" .to be "meaningful-" 

Succeeding sections of this report will- present many of the_bas1c 
principles and. techniques:. involved .in thermal math- modeling^. Experience, 
ofc-course,-can. only, be -acquired- from hands -oo_fam11.iartty with-feal themal 
systems and. participation Jn the- modeling, and analysis thereof-. Engineering 
Judgement can- probably-.be described more. accurately as -the. result of ab- 
stracting from- the-body of unique familiar information, a general under- 
standing which can be. extended to guide the -investigation and - comprehension 
of new. and unfamiliar areas. As such. Engineering Judgement cannot be 
presented in a table, ora figure, or even an entire-book. However, it is 
possible to present a collection and discussion of examples of thermal math 
modeling which will serve, not as a source of practical experience, but as 
a body of unique and familiar information from which the reader may abstract 
as_much Engineering Judgement as he Is able. Such a collection Is included 
in-the latter half of this report. 


2".0 THERMAL MATHEMATICAL MODELING 


2.1- Network Solution 

Two systems- are said to be analogous when. they both have similar- equa- 
tions and boundary conditions; and the. equations describing, tiiebehavlor of~ 
one system, can be transformed 1 IntCLthe equatl ons_ f or ths-oth ir by simply 
changing .symbol* o£ the_var1ables. Thermal and erVectrioal j.ys terns are two 
such analogous systemsas shown- In table 2-1. 


tjuantlty 

Thermal 

Electrical 


System 

System 

Rotentl al 

T 

E 

Flow 

$ 

I 

Resistance 

R 

K 

Conductance 

G 

T~ 

R 

Capacitance 

G 

<r 

0hm 4 's Law- 

Q_= GT 

,—T 
1 IT 


Table 2-1. Thermal-Electrical Systen Analogy 


The analogy between tbermal_and electrical! systems' allows the engineer 

to-uti I tae th e wl dely known basic laws-such as Ohn^s Law and Klrchhoff^s 

Laws used- for balancing networks * numerical techniques used to solve the 
partial differential equations describing such •systems have been conveniently 
adapted to computer solution*, thus-enabling thf* engineer te readily compute 
temperature distributions and gradients- ef complex-physical thermal networks. 

Thermal" analyzer computer programs-have-bien developed- which require 
the-user to define a- thermal— network- of" tho s^stent-anaTogous to an elecr- 
trfcal circuit. The network components - fire fiput Into the computer and 
pre-programmed routines perform the transienl or steady state-solutions." 











This secttor»_di s cusses the development-of- a thermal-network .and the 
numerical techniques for solving, this network. 

~2 .2. ModeHng_Elements 

2.2JL Nodes 

2; 2.1 .1 Concepts 

In order to-develop.a thermal network and apply numeti cal' .technl ques~ 
to- its salotfon % .iie-1s-necesssiy to. subdivide the thermal" system Into a: 
number of finlte subVolumes- cailed-nodes. — The thermal' properties-of" each- 
node are constdeced to Jjeconcentrated_afr-the_ceritral nodal point of- each-, 
subvolume . Each. node_ represents twa_thermar network elements »_ a tempera* 
ture (potential) and a capacitance (t hemal mass) as shownJn. Figure 2.1 , . 


hti Cno<M i) 
T 2 fnOde 2) 


^ W 

Tj.Cs {node 3) 


Figure 2-1. Nodalfesatlon >s?T - 

The temperature* Jr,_ass tgned- to a uode »epresents-the average mass - 
temperature-of the sub volume. The capacitance* C* ass1gned--to a node fs - 
computed f rom-fehe thermophysi cal- properties o t the subvolomeHnaterialT 
evaluated at the .temperature of^the node and Is assumed-tb be -concentrated 
at the nodal center^-ofthe sub volume. Because a-node "represents a- lumping 
Jjf parameters t€t«inol« point in -space* the temperature distribution 
thlough-the subve1uroe-4mplied"-by- the nodal— temperature- is ttnew-as shown- 
In-Figure 2. 2e-and-not a step function as illustrateddn- Figure^ZUh, — — 
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(a) 

figure 2.2. 

In a homogeneous material, the temperature at a point other than the 
nodal point may be approximated by interpolation between adjacent nodal 
points where the- temperatures are known. 

The error Introduced by dividing, a system into finite size nodes rather 
than vo l ume s dx^ where_dx approaches zero is dependent on numerous consid- 
erations: Materi aL thermal properties^ boundary-conditions* node size, _ 
node center placement, and time increment-of transient calculations. The_ 
techniques for- proper nodallzation-to -minimize- the error will be discussed' 
ln-a-3 ater section... 

2. 2.1,2. Node Types 

To this poinfonly-nodes^whlch represent subvolumes with a finite 
thermal mass {capacitance) have been discussed. In many instances-, two 
other types of nodes are required to define a thermal network. They are 
nodes having a zero- capacitance or an-4n£in1 te capacitance. Thermal' 
analyzers such as SINDA usually- gitfe the three types of nodes particular 
names as follows: 


(b) (c) 

Temperature Distributions 



Name 


Node Type 

Einite Capacitance- Diffusion 

Zero. Capacitance Arithmetic 

Infinite Capacitance Boundary 


The. diffusion node has a finite capacitance and Is.. used.. to repre= 

sent normal material. nodal Izatloru It Is- characterized by a gain or... 

loss-of potential-energy vdrl eh .depends on the capacitance value-, the net 
heat flow into the node, and the time over which the. heat fs-flowlng. 
Mathematically, a. .diffusion- node is defined- by the expression 


The arithmetic node (zero, capacitance)' Is a physically unreal quantity, 
however, its effective use with numerical solutions can often be helpful In 
Interpreting results tii such applications as surface temperatures, bondline 
temperatures, and node coupling temperatures. It also finds use In repre- 
senting thermal system elements which have- small capacitance values In 
comparison to the large majority of the other nodes *n the system which 
results In computer run time reduction with minor changes in overall 
accuracy — Examples of these could be small components such as bolts, 
films, or-fillets; gaseous .contents of small ducts or-tubes; and-low mass 
Insulations^ Arithmetic nodes should be few In number when contrasted to 
the total number of nodes ia.-the_network* _Ihe temperature_of an arithmetic, 
node responds Instantaneous 1y_ta_its surroundings. Mathemat ically, an 
arlthmetic node Is defined by the expression; 

EQ = 0 

The boundary node (Inflnite-capacltanee) Is used- to represent constant 
tempera ture^sourees within a thermal network - . Common uses are t_ deep 
Space sink temperature-, recovery temperature - , lunar surfaee temperature, 
in addition* a-boundery- node may-be- used te represent thermal system com- 
ponents such as the bulk propellant in a Targe” tank whieiv has -a very- 
large* thermal mass (capacHaneej- Mathematically a- boundary node Is 
defined as: 

T # Sons teat 


6 


2. 2. 1.51 Method of- No dal 1 zation 


The placement of dlffyslon node centers and-the.cho1ce_of node.shapes 
Is dependent on several factors: (1) the points where temperatures are 

desired, (2) the- ejected temperature distribution., (3) physical- reasonable- 
ness',. and (ALtbeeaseof computation. The- actual size .of- the node Is 
dependent on other considerations: 0) accuracy desired, (2) structural 

design. (31. computer storage-capabllttlesr, and- (A-)- computer time required. 
Eachi-factor, however^ embodies other-eonsiderati ons. For example, to 
anticipate the.expected temperature .distribution one must draw heavily on 

En gine ering Judgement: as to the effects of “the-expected boundary. conditions. 

and associated material properties. 

In general* the shape of a diffusion-node is chosen to be- a simple 
geometric figure whose areas and volume can be easily calculated. Irregular- 
ly shaped structural members may- be approximated by_simple shapes by em- 
ploying assumptions that are consistent with the desired results. Mode 
centers are assumed to be located at the.mass centroids of the nodes. In 
some cases, nodal divisions are decided first, with the node center locations 
being thusly defined as a consequence. In these cases, nodal -boundaries 
will usually lie along structural boundaries, and- structural members will be 
divided in_a symmetric and. equal-fashion* In other eases, output require- 
ments will dictate the locations of node centers, with the nodal boundaries 

assigned^as a consequence. These two approaches are IllustrateiUn Figure 2-3* 

In the- first case. Figure 2-3Ui.1t ts-deslred-to prepare a general model 
offhe-s t-ructure, but In the second case, Figure- 2^-3(h) it is desired to- 
modal the response of the two thermocouples located on tee bend line between 
the two members. 






(a)Cftriaral Node-Boundaries 


(t>) Node- Banters at thermocouple Locations. 


FI gure 2-3v A1 temate Nodalizati on JSelhods 

The above example. aUudes to a general desirability for. rectangularly 
shaped nodes* This is true. for the_slmple reason.that.it is easy to. com= 
pute the areas and volumes required for the Input calculations- Such 
simple nodal shapes are in keeping with current engineering practice. 

By contra st, Dus tnberre (reference 1) suggested that. nodal tzatl on be per- 
formed in such a manner that the-paths of heat flow assumed a triangular 
pattern, as shown in Figure 2-4(a). The only drawback to this theoreti- 
cally sound approach Is that it remafns for the engineer to compute the 
volumes of the irregular polygonal nodes which are the consequence of such 
a tact, as shown in Figure 2*4(b). 


(a} Trtangulap Heat 
Flow Paths 
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fc) SlmpJa Rectanqular 
NodaUaatlon- 


Figure 2-4. Polygonal Noda1i2ation Vs Rectangular Nodatlzation 

Note how much simpler Is the rectangular nodallzatten approach indicated! in 
Figure 2-4(e). As might be expected# to- achieve- the Sana simplicity' of 
calculation# circular structures are- nodall-zod In pie-shaped Wedges#- Annu- 
lar rings, or a combination of the- two, as shown in Figure 2-5. 



I ) 
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e 



Figure-2- 5. Nodalization of Gi rcular- Elements— 


Boundary nodes_are^used. to define points J toes , or surfaces ol 
constant temperature fn one* two, or- three, dimensional. models respec= — 
tlvely. The-physi cal' location of~a boundary node Is- determined solely- by 
the conduction paths connected to Hr A. single boundary node may be used 
to model all. boundaries which are at the same temperature. Thlsr point Is 
Illustrated -In Figure 2-S which shows, that the Indicated boundary node- will 
suffice as a model of the entire* constant temperature edge. of the structure 
(In this ease, at 100'T). 



Figure _2-iu Sample Boundary node 

Arithmeti erodes' have-a numberrof uses which a nr consequences of, the- 
fact-that-such: nodes^serve i^-an engineering -model the proverbial "wafer - 
of thlckness^dxvJwhepe- djc-approaches zero".— A- typical application lies- 
in the-rodellngof^exterior surfaces^of^reentry vehicles which are— 
oftsn-subjeeted— to severe, rapidly changing,- boundary conditions, -to the 
physlpansystem>_the-surface-temperature remains very dose- to radl-at ion- 
equilibrium with' the surface^-heatl ng-ratModlcatl ng that- -thie-sy st em -ean- be 
accurately-simulated by the-use of^-a -surface arithmetic node. This appli- 
cation H Illustrated^ Figure- 2-T. 



Figure 2-7. Use of Arithmetic Nodes— to Model Surfaces 

The case where heat flows- from a surface by conduction is usually one 
in which two structures are bonded together and a bondline temperature 1s- 
s ought. When the- structures are homogeneous ,-a bondline temperature may be - 
established by simple linear interpolation between. the nearest node centers. 
When the materials are dissimilar, It Is: more appropriate to use an arith- 
metic node at the bondline, leaving to the compute/ the process of per- 
forming a conductance- weighted-average of the ed jol ntng_dlffusi on node, 
temperature which* In- essence, 1s_the res ultroffl nding-the steady state 
(heat tn. * heat out ) tem perature for an arithmetic node. 

Arithmetic nodes may also be. used advantageouslyrln-piace-of diffusion 
nodes-Whieb have a capacitance- that is- small when- compared to the great 
majority of nodes in- the system* This -often -occurs when modeling a small - 
quantity of-gas tn a tube or other enclosure, or^when modeling small struc- 
tural parts, such - as>4res, bolts, fillets, films, and' sheets-, where 
detailed temperatures are desired (which precludes lumping sueh-ttems along 
with— larger nearby nodes) - . The correet use of arithmetic-nodes in these 
cases generally results in a considerable saving- of computer time when the 
model Is processed. 



2.2.i.4 Computational Methods - Nodes_ 


In developing. a thermal network, computations with respect to nodes — 
are generally limited to- calculating the capacitance of diffusion nodes*. 

The following formula is used* 

C- p • V - Cp. 

where: 

C = thermal capacitance. BTU/^f 
p ■ density — LB/ FT 3 
V » ... volume-^ FT 3 

C p _=- specific heat - BTU/tB-^E 

The specific heat (Cp)„ and the density (p) of materials may vary with 
temperature. The necessity to utilize temperature dependent-properties for 
an analysis depends on the degree w1th~wh1ch_the properties vary and the 
temperature range over which the capacitance of- the material will be calcu- 
lated. Most thermal analyzers can accommodate temperature varying thermal 
properties. 

The use of arithmetic nodes may also require some computations _ Re- 
placeraent-of small capacitance diffusion ncdes with an arithmetic node must- 
be preceded by computations to. verify that- the capacttance-conductor effects 
are such tbatulhe node-in question wfll—essentlally reach steady state 
temperatures during the_ttme step requited by. the larger nodes. The-use of 
an arithmetic-node -to predict surface temperatures, where surface-radiation: 
ocu/ery bigh Jjeatingj»tes-are 1 nvo 1 ve d'requ ires.“xa ref ul -analysis to insure— 
the- stability of- the arithmetic^ node. — Stability- eri teri a and solution — 
techniques - are- discussed 1 n-section 2 . 3. Prom^this section* It: can be seen 
that _ solution-teehniques using: linearized ^last-pass" temperature Values* 
may requirethe use of* analyzer- control- constants to restrict"the maximum 
node temperature change on- computation- time step. The-engineer must further 
bo cautioned against using, coupled, arithmetic nodes- without -a complete 
understanding of the impltcatlons and requi red-analyzer-control constants 
used-to insure a vallCTsolutl on. 


t. 2.2 Conductors 
2. 2. 2.1- Concepts 

Conduetors-are-Jthe- thermal matb_modeUjig_ne.twork elements which, are 
used to represent the heat flow paths -through which, energy is transferred 
from-one-node-to- another node, figure 2-8 Illustrates the element-node 
temperatures CrWcapacitances. (C) and conductors U) which comprise a 
thermal networks — 


8 T 

V W • 

r 

Figure 2-8. Thermal Network-£lements_ 

2. 2. 2. 2 Types of Conductors 

The three- processes by which heatuflows_from a-reglorrof higher 
temperature to -a regton-of lower temperature are-conduction, .convection'* 
and radiation^ Conduction Is -the. process by whim-heat flows w1thin_a 
medium-or between different-mediums _in -direct-physical contact. The energy 
is' transmitted by molecular-communf cation, figure £-9 1 ULustrates_ con- 
duction conductors. _ 




7 

- 

r 

i 


rJiA/L 

- 




— ” 




- -- 


— Figure 2-U.- Conduction-Conductor" 

Convection is the process of energy_Jtransport by combined^ctidn of heat 

conduct toiv-energy storage* and : mixing, motl on> Meat will flow by conduction 

from' a-surfaceto adjacent particles of fluid;— thett-the^-flufd particles Will ) 
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move to a region of lower temperature where they will- mix with* anti transfer 
a part of their energy te, other fluid particles. The energy Is actually 
stored in, the fluid particle! and Js carried as a resu'.t of their m ss 
motion. Eigure 2-10. illustrates the-conveetlon conductor* — 



Figure 2-10. Convection Conductor 

Conductors which represent conduction or convection paths are re- 
ferred to as linear conductors because the heat flow rate is a function of 
the temperature difference between nodal temperatures to the first power. 

' !' Q = 6 (Ti - Tj) 

Radiation Is the process by which heat flows between two bodies when 
the bodies are separated in space. Energy Is transferred through electron- 
magnetic wave phenomena. Radiation conductors -*re- termed non-linear be- 
cause the heat flow between two- surfaces by. radiation Is a function of the 
difference of the fourth powers of the surface temperatures . 


0-6 (Tt 4 - - Tj 4 ) 


Figure 2-11 Illustrates the- radiation conductor. 
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Figure 2.-11. Radiation Gonductor 

Fluid flow thermal systems may also be simulated by thermal modeling. 
Energy stored In the thermal nfess (capacitance) of a fluid lump (node) Is 
transferred from one point to another by the movement of the fluid mass... 
This type* of conductor Is generally referred to as a mass flow conductor 
and Is Illustrated In Figure 2-12. The mass flow conductor Is also 
linear 


Q - 6 (Ti - Tj) 



FI gure 2-12. Mass Flow Conductors 



2.2.2. 3 Computational Methods - Conduction Conductors 


Conduction conductors are confuted from the equation: 


where r 


_ kA 
s - CL 


G j*_thermal conductance - BTU/HR--E 
k = thermal conductivity - BTU/HRtFT-T 
A *- cross sectional area through whleh heat.- flows* * FT 2 
L = lengtiLiffitwsen .adMnlng nodes ~ FT 

The thermal conductivity (k) Of materials may. vary with temperature 
or other influencing factors within the system* the cross-sectional area 
through which the heat flows- (A) and length between node centers (l) are 
determlned-by the size and shape of the adjolnirtginodes. As with the 
capacitance calculations, the necessity te use temperature dependent- 
properties depends on the degree with whteh the conductivity changes over 
the temperature range expected during the_analy$1r,. 

2t 2.2.3.1 Rectangular Nodes 

the length* L, of the heat flow path*, used for conduction conductance 
calculations for rectangular nodes Is Hie distance be tween-node- centers, 
and tiie area* A, to be used- Is the: cross-sectional area. perpendicular td-the 
line joining the node centers. The convention- Is depicted in figure 2-13, 


Figure 2-13. Simple Conductor- Representing a Heat 
Flow Path through Material 


( 





2.2.2. 3.2 Circular Sections 

For cortductOrs~between nodes Which are circular-sections* the con- 
ventions shown In Figure 2-14 should be used: 


A _ Qd 

C TnTrJ7rfT 

*here: 

e - radians 
d - FT 

r - consistent units with r 4 
o i 

r-j - consistent units with r 0 



Figure 2-14. Area and Length Equivalents for 
Circular Section Nodes 


2.2.2. 3.3 Parallel Conductors 

Two or -mere- parallel conduction paths between nodes may be summed to 
create- one- conductor -Valuft_by the fol 1 owl ng- equation r 

G-pa-Gj + 62 + 6^ 

This may: be~helpful in compyt1ng_an_equ‘Ualent conductor between two nodes 
as Illustrated in Figure 2.15. 





h 


s T * 6? + 62 

Figure' Z. 15, Parallel Conductor Flow Paths 


2. 2. 2. 3. 4 Series Conductors 



Two or- more series, conduction paths between nodes may be combined to 
create one conductor value by the following equation: 


1 


i 



* • * 


Gfj® 


T - — 1 



This may be helpful-in compUti ng_ the conductors-between two dissimilar 
shaped or dissimilar material nodes~as shown in F4gure_.2.16_ 




G i - 


k i A i 



1 - 


®T " 1 TT 


G 1- G 2 

spnq 


Figure 2. IfT. Series Conductor Paths 
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2 2.ZA Computational Methods "—Convection. Conductors 

Con*?ect1on conductors are computed from the expressions 


hA 


where: 


G = thermaV-conductance_=^ BtU/HR— °F~ 

h = convecti ve-hea±transfer_coe£fi ci'ent_r BTU#I&-FT 2 - <> F~ 
A>-.surface .area in -contact with the fluid. - Ft 2 

&_is~the -product'of the average unl trthennal convecti ve conductance 
h {convective_heat transfer or film coeffi cienfl and -the nodal surface area _ 
A in contact with the fluid. However, h Is a compl fcated^unetion of 
fluid flow, the thermal properties or the fluid medium, and the geometry 
of the system. 

Since the convective process of heat transfer is so closely linked to 
fluid motion, it is first required to establish whether the fluid flow is 

laminar or turbulent. In laminar- flow the fluid moves in layers and the 

fluid particles follow a smooth- and continuous path. — Heat ts transferred 
enly by_moleculad-eonduct1on within the fluid as well as. at the interface 
between the- fluid and the surface. In. .turbuient-fTow the path of the fluid 
particles is' Irregulars- and although-the general trend-Cf the motion tSLln 
one direction*, eddies or mixing currents -exists In addttlon-to the con- 
duction- mechardsoLbeing modi-fied, incre a$e d_h eai-iransf er occurs in turbu- 
lent flow when -energy Is carried- by fluld-pariUles-across flew streamlines 
and mixes wfth-other parades of the^fluid. 

In addl tlon. Jtodcnowl ng- whether the -fluid motion U laminar or turbu- 
lent", It -+s necessary to know^ the prOcess-by which- the motion was Induced. 
When-4;He heat flows between- the fluid and-fhe surface as a result of fluid 
motion -eaused by differences- In fluid density-resulting from teir^erature 
gradients in the fluid* the heat transfer mechanfsBrls called free or 
natural- convection. When the motion is caused by some external agency*~such 
as a pump^-Ot-blower , the heat transfer-mechanism is cal led-force d^onvention. 




■CD 

The following table Illustrates typical-values of average heat 
transfer coefficients encourttered-tn-jengi neeri ng practice. 


Convective Medium 

Convective Heat-Transfer Coefficient- 
h - BTU/HR-FT 2 --°F 

Air, free. Convection. 

1-5- 

Air,. Forced Convection 

5^50 

Oil*. Forced Convection 

10-30CL 

Water, Forced Convection 

50-2000 

Water, Jotllng. 

500-10,000 

- Steam, Condensing.. 

1000-20,000 _ 


Table- 2-2. Order of Magnitude or Convective Heat Transfer Coefficients 

Equations for the computation Of conVectl ve. heat transfer coefficients 
will be divided into three categorlesr natural conveetion, forced convection 
in tubes or ducts, and forced convection over exterior surfaces. The equa- 
tions presented for the calculation of convective film coefficients are 
the most generally used expressions. Others are available and are applicable 
for many specific thermal problems- or analyse,. It should be remembered 
that the predicted values for h are only approximate. The accura<y of the 
heat- transfer coefficient calculated from any available equation or graph 
may be no better than 30 percent. 

2. 2. 2. 4.1 Combined Natural and Fereed Convection 

ljt_Cases where both natural and forced convection are combi ned+_i tils 
left te the engineer to decide which- heat-transfer phenomena Is significant 
and utilize tee proper equations to compute the- effective film coefficient* 
(h). The -following deternvtoations^hould be helpfuTr 

1) Calculate the Grashef number (Gr) 

p 2 g MT - T,)'t3 
Gr 



where: 



Gr ■■ Grashof' number- - dimensionless 

g —gravity constant - 4.17x10^ FT/HR 2 

0 “ coefficient of volumetric. expansion 
@ film-temperature-- 1/°R 

T> temperature of the body, surface — °F 

Too ■- temperature of the surrounding medium-°F 

L“ length -of flow p_ath — 0“ 

a* fluid Ad scosity. — LU/HR-fT 

2) Calculate the .Reynolds- number -(Re) 



where: 

Re * Reynolds number - dimensionless 
U ■ fluid velocity - FT/HR 
p = fluid density - 10/FT* 

L * length of flow path - FT 
v * fluid viscosity - LB7HR-H . 

3) If Gc_< 0,225 Re 2 the ef.fectiof natural convection on the_ 

avecage heat transfer coefficient-for 
pure-forced convection fs less than 5-% — 

4) if- Gr > -KLO Re 2 forced- convection has. negligible- effect 

on-natural* convection— 

5) In the region where both, free and forced convection effects are 
of 1 tile same order-of magnitude, heat transfer is increased by 
buoyancy effects acting in the direction of-flow and-decredsed - 
when acting In the opposite direction _ 

6) In cases where it is. doubtful whether- forced Or free convection 
flow applies, the heat- transfer eoeffictent-ls generally calcu- 
lated by using forced and- free- convection relations separately 
and the larger one-ls used. T-ho-aeCUracy of'this ruin of thifinb 
Is esttmatedrto be about-2535. 
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2.2.2.472. Natural? Convection Equations 

Natural convection occursr when fluid Is set In motion as a result of 
density- differences due to temperature variations in the fluid. A great 
deal of ■ materi al has, been published on heat transfer coefficients^ result 
ing in many, techniques, with a wide. variety of results., leaving, the 

"occasional user^nof the= Information easily .confused. Presented herein 

is- a wide variety of applications condensed i nto a_smal 1 nunber of groups 
The 1 nl tial data. was. obtal ned from the Oil and Gas: Journals Equations _ 

used are. the ones which. appear to_have the wl dest acceptance .for-the 

particular selnof cond1tlons..involved 

Correlations of natural convection heat transfer usually take...- 
the formr 


h = K (Br Pr) n 

where: 

h = heat transfer coefficient 
K « an empirical constant 
Pr * Prandtl number 
6r a Grashof number 
n * an empirical exponent — 


to apply the above equations-to an actual problem solution, the 
following steps should be performed: 

1) Calculate the product of -th* Prandtl and Grashof numbers 

p 2 e C p g aT b 3 
6r-Pr » 


where: 


u k 
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k * thermal- conductivity 

0 film temperature ~ BTU/FT 2 -HR-°F-FT 

p « (tensity 0 film temperature - LB/ FT 3 

& = coefficient of volumetric expansion. 

0 film-temperature — l/^R 

g * gravity constant ^ 4.17*10 8 FT/HR 2 

C p «*■ speclflc-heat .0 -film, temperature - BTU/LB-°F 

y viscosity Q-film temperature - LB/FT-HR 

L » h eight "or length of surface - FI 

AT - temperature difference between the wal Land and -fluid 

Film temperature » Ttfa ^ * T flMld _ op 

C t " 

2) If 10 3 < fir Pr < 10 9 assume flow Is laminar and go to Step 3. 
If Gr Pr > IQ 9 assume flow Is turhulent and go to Step 4. 


3) Laminar Flow 

h » C] | (fir Pr) 1/,f 

where: 

h “-heat-transfer^coefflcient - BTU/FT 2 -HR^-°F 

k * thermal conductivity atifilm temperature- - BTU/FT-HR-^F 

U- height or length of surface - FT 

Gr Pr * Product of- Grashof and Prandtl numbers 
from above - dimensionless 

Cl is dependent on the geometry of the system and can 
be determined- from- Ptguce^-2-17 

• T wan + ^fluld 
T ^ 


Film temperature 


4} Turbulent Row 



h » Ci_ £ (Gr-Br) 


1/3 


where: 

h » hea transfer coefficient ~RU/fT^JiR*±£ 

k * thermal conductivity at: fiJm temperature BTU/R'-HR-^F 

L *-hetghtrojr-iength ofisoirface — R 

Gr Pr * Productaof Grashofand Erandtl numb ers 
from- above - dimensionless 

Cv H dependent on the geometry ef_the system and can 
be_detera1ned from-F1gurd-2-T8 


Rim temperature « 


Wl * T floid _ 



o 
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laminar flow 


® $ ^ Wara Sas tell ^ p u *& 

^ N 

Short tell <2 ft -Short Wall < 0.15 ft. 

Cl " (L548 c t *. 0;54& 


cool 

Wall Liquid 


Short Wall <0.15 ft 
Ci- = 0.55 


Cool Gas 

t. j 



Warm 

Wall 


Short Wall. < 2 ft... 
Ci " 0,548 


Cool Wall 


Small Wall < 2 ft sect Small Wall < 2 ft sect 
Cl ” 0,35 Cl a 0.71 


Warm Wall 


r 

Cool Gas 


Small Wall < 2 ft sect ( Small Walt < 2- ffc_sect 


Ci ■ 0.7-1 


G'l * 0-36 


Cool Tube 


Small Tube: D < l_in 
c l-~ 0>53 Sub D for L 


f— Warm Gas 


Cool Tube' 


Gas or i 
Liquid >JF 

cf 

& 

f Cool 


— f 

C Warm 
Liquid 


)L_ Liquw ^ Warm ^ 

StelT Tub* lln < B < 4 1„ MoT *»»• 

Cl ■ 0.47 Sub D far I ,, , s ™^' Dtmetef 

L Liquid: G< 0.15 ft 

Gas: R < o.50 ft 

Cl - M3 Sub K for L 

Nowre 2-17. E-pIftM Con«t.nf, foHi.tu«l Cb^ctloh - L«,n.r FTtm 










TURBULENT FLOW 


lool \ 
tell •. 



Warm Gas 


Long Wall > 2zft 
C, * 0.1? 


Cool \ 
Wa11\ 


Warm 

Liquid 


Long Hall > 0.15 ft‘ 
C t *_ 0,13 



Cool Gas 


WarmV^- Cool 
Wall \ WLtquI 


Long Wall_>-2 ft ' Long-Wall > 0.15 f 
Ci - 0.13 Ci^O.13 


Warm Gas_ 


Cool Wall 


Warm-Liquid 


Cool: Gas: 


Cool Wall 

Large Wall > 3' ft 2 


Warm Wal 1 


Gool 8as 

Lange WaVl > 3 ft 2 


Warm Gas _| Cool WalL. Warm Wall 

Large Wall > 3 ft 2 Large Wall > 0.25 ft 2 Large W&U-> 3 ft 2 

62 Ci « 0.08 


Cool Liquid I ^ 


- / ///// 

Warn Wall 

Large WalV> 3 ft 2 

Cj » 0.162 


^ r 

Warm Liquid 
Cj « 0.128 


Coot Wall 


Cool Liquid 


j^^jcooi Tube 
Warm Liquid 

Ci - 0.128 


Warm Gas 
Cool Tubety 


Large Tube: 0 > 6 In 
Ci - 0.11 Sub D for L 


Warm- Tube 


Cool-Gas 
Ct - 0.128 


M.. . . ... Warm or Gas or 

| warm Liquid Cool sphere * Liquid 


Cool Tube' 1 


Large Tube*- tr> 6 In 
Ci « 0.11 Sub D for L 


L-arge Sphere f 

Liquid-: R > 6.15 ft 
2811 R > 0.-50 ft 


Ci - 0.15- Sub R for t 

Figure 2-18. Empirical Constants for Natural Convection ~ Turbulent Flow 














2. 2. 2. 4. 3 Forced Convection In Tubes and Ducts 




TO-a pptv the. empirical relationships For forced, convection heat 
transfer eoeff 1 clents in tubes^and ducts, the fa! lowing .steps should be 

performed. Comments on thermophyslcal property, determination and 

entrance effects* are Included at the end of the section. 

1) Calculate the hydraulic diameter. Oh, which is defined as: 


. flow cross-sectional area 
°H * 4 wetted perimeter 

the hydraulic diameter for some common shaped tubes or ducts is: 

Round: Da m D» where D is the diameter of the tube 

Square: Dh ■ b where L Is the length of a side 

Rectangle: Du a » where a and b aretha lengths 

a+t) of the sides 


Elliptical: Dh 


4 


abVa 2 + b 2 


where a and b are major and 
minor axes 


Annular: Du * Do - Oi,- where D 0 Is tbe-outeit tube diameter 

Df is the inner tube diameter 

2) Compute tht^ Reynolds number. Re-, from the^fol lowing equation: 


He r Oils. 

where: 

Dh » hydraulic diameter, defined above - FT 
U ■ velocity - FT/HR 
p ■ density-- LB/FT 3 
u » viscosity - feB/FT-HR 

3)- If Re < 2100 the flow can be : considered-laminar, go to step- 4 
If Re > WMJOO- the-flow-cah- be considered turbulent* go- to Step S 

If 2100* < Re > 10,000 the flow can be considered transitional, 
go- to step 6- 


) 


) 
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4) Lanvlnay* FTnw 


d 




1 .86 k ( Re P* — \ 

V Dh 2 U 


1/3 


fi = 3.66g is a minimum for 1 amfnar flow in long-tubes; 
for very short tubes 


(fejf «-00S R«i r- 


Re Pr k 
“ 5C 


In 


t- 


2.654 


167 (Re Pr Dh/l) 


• 5 ) 


where : 

h » convective heat transfer- coefficient - BTU/FT 2 "HR-°f 
Re- » Reynolds number from step- 2 - dimensionless 
Oh * hydraulic- diameter - from step 1 - FT 
L_»r_ tubfr or duct lengths Ff 
B = tube, or duct- diameter - FTT 

C 

Pr - PrandtT number » - dimensionless 

C p ■ specific heat - BTU/LB-°t= — 
k * thermal conductivity - BTU/FT=HR~?F 
w-* viscosity - LB/fT-HR- 
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5) Turbulent Flow 


k n *® r, * 33 

h - 0.023 ^ Re Pr 

Wher8, h « heat transfet^coefflclent-^ BTU/FT 2 -HR-°F 
k * thermal conductivity.- BTU/FT^HR-°F 
D h ^hydraulic-diameter- from step! - FT 
Re » Reynolds- number - from step 2 dimensionless 

Cp v 

Pr-= Prandtl -number ° ^ — - dimensionless 

Cp. e-specific best .- BTU/LB-°F 
y * viscosity - LB/Ft-Hk 

6) Transition Flow i/3 t v b \* 1H 

h ° Ci Cp U p Er ^ — I 

where: 

h ■ heat transfer- coefficient - BTU/FT 2 -HR-°F 
Cp * specific heat - BTU/LB-°F 
U *» velocity - FT/HR 
p «- density - LB/ FT 3 

C ti 

Pr = Prandtl number * -jr— ~ dimensionless 

u = viscosity evaluated at the average, film 
temperature — LB/FT-HR 
uK- viscosity, evaluated. at the built fluid 
temperature - LB/FT-NR 

uc = viscosity evaluated at the s urface temperature 
LB/FT-HR 

R * thermal conductivity - BTU/FT-HR-*F - 
Ci * a function of Re and can be- determined- from 
Figure 2*19 
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FORCED CONVECTION TRANSITION FLOW IN TOBES AND DUCTS 



Thentiouhvsl cal Properties 

Because_the film coefficient (h) Is a function of the themophyslcal 
properties of the fluid, It. is Important to account fo»- any significant 
variation In the properties as a result of the temperature- of the fluid. For 
liquids ^generally only the temperature dependenc-e of the viscosity (v) Is of 
major Importance. For gases, the other properties density (p). conductivity 
(k) and specific heat (G p ) usually vary significantly. It is recommended that 
all of the properties except" Cp be evaluated at the average film temperature 

of the fluid Tf defined as: 

T f * 0.5 (t s + Tfc) where T s ■ surface temperature 

Tp « fluid bulk temperature 

C p should be evaluated at the fluid bulk temperature. 


I ) 

All of the above equations, with the exception of the one foe com-v 
puting h for laminar flow tn very short tubes, ignore antranee effects on 
the film coefficient. Thr Importance of entrance effects depends on thfr 
fluid flow condition and the length (L)/hydraul1c diameter (Dh) of the 

tube or duct. Entrance effects are appreel able where: 

For laminar flow ^ < 50 

FOr turbulent .flow t- < 10 - 

H 

In general, entrance effects Increase the effective film coefficient. 

The local heat transfer coefficient (h x ) divided by the free stream 
heat transfer coefficient (h») approaches 1 as ^ Increases. If the detail 
of an analysis requires the consideration of entrance- effects on the con- 
vection coefficient, literary research may be helpful In finding an equa- 
tion approximating the effect. . An example of entrance effect correction 
is as follows: consider turbulent flow in short circular tubes where 
(2 < L/D < 60), the effective film coefficient can be approximated by: 

.7 

h E = h + h (£) where 2 < L/O < 20 

h E ■ h + h (1^-) where 20 5 L/D < 60 

where:* 

h> - effective film coefficient corrected for 
E entrance effects - same unit's as h 

h » turbulent film coefficient calculated from step $ 

\ 

| D - tube diameter - units consistent with L 

L = tube length - units consistent with & 

1 

; { 
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2. 2. 2. 4. 4 Farced Convection Over Flat Plates. 


Emperlcal relationships for forced convection ever the exterior surface 
of a flat plate are given_fbr bath- the entire length (14 of the_plate -and at 
any .intermediate length (x). Equations for both 'locations are given fn step=_ 
wise fashion r 

1) Calculate -the Reynolds number Re 




where: 

U» » free stream velocity - FT/MR 
p = fluid density - LB/FI 3 
vi = fluid viscosity - LB/FT -Hit 
x ■ Intermediate plate length - FT 
L • total plate length - FT— - 

2) Tf Re < 5xlQ- 5 flow. Is laminar* go to-step 3 
If - Re > 5xlC 5 flow is turbulen tv-go to step 4 

3) Laminar -heat transfer coefficient" (h) - 

k 1/2 . 1/3 

Evaluated at x* h = .332 ~ Re x Pr 
Average, val ue for plate with length L» 

k 1/2 i/ 3 

h » .664 £ Re L Pr-' 

where: 

h » convective heat transfer coefficient - BTU/HR fc .FT 2 -°F 
k * thermal condUotivi ty - BTU/HR=F?-°F 


X - intermediate p» ate length - FT 
L * total plate length .—FT 

Re ^.Reynolds' numbe r front-step 1 ~ dimensionless:- 

C« u 

Pr * Prandtl number a dimensionless 

y * vlscosi ty -LB/JrT-HR 
4 L Turbulent-heat transfet coefficient (ft) 

k n * ® 

Evaluated at h "■ 0»0288 — Rex p r 
Average value for plate with 1 eng tit L, 

t * 8 « l / 3 

h “ 0.036 £ Re L Pt_ 

where: 

h - convective heat transfer coefficient - BTU/HR-FT 2 - 

k * thermal conductivity - BTU/NR-FT-°F 

x “intermediate plate length - FT 

L o total pi ate length — FT 

Pe «_ Reynolds number from step -1 -• dimensionless 

C y 

p r « Prandtl number * 7 “ ~ dimensionless 

y ■ viscosity - UB/FT-HR 


g. 2, 2. 4.5 Forced Convection Oyer,, Cylinders 

Empirical relatlonsh-lps^-for-flow over cylinders are given for the 
stagnation point* -forward-portion, of the cylinder, and for laminar flow over 
tiie total “cylinder. 



Figure 2-20. Flow-Over Cylinders 



Stagnation Point 
wherer 

h « ecnveetive heat transfer coefficient - BTO/TT^-HIU^F 
k-« thermal conductivity-- BTU/FT-HR-°F 
I]* * free stream velocity - FT/HR 
D «* Cylinder diameter - F7- 
P » density - LB/FT 3 
p » viscosity - 1-6/ FT- HR 

e n v 

Cl - Is a function of - Pir * -fr— arid can Be approxl mated 
from Figure 2-21 

Cp «• specific heai-- BTU/tB=*F- 

33 


Figure 2-21. Empirical Constant for Cylinder 
Stagnation Point 


Forward portion Of cylinder where (0 o< 8 <80°) (see- Figure 2-20) 



where: — 

h * convecti ve heat transfer -coefficient - BTU/FT^-HR-^F 
It * thermal conductivity * 1JTU/FT-HR-°F 
D * cylinder diameter-- FT 
tf «5 * free stream velocity - FTyHR 
p * density - LB/FT- 3 
y * viscosity ^-LB/FT-HR- 
Pr «■ Prandtl number «— |p- * dl mens io«tl ess 
Cp « specific Heat - BTO/tB-°F 

p » angle measured from the stagnation point - degrees 


Average fil.n coefficient for laminar flow over a cylinder 



for a gas h = C | ( fe P P . ) 


n .31 


for a liquid h * 1.1 C p r 


h convective heat transfer -xoeffl c1enfc-BTll/FT 2 -HR-°F 
k * thermal conductivity - BTU/fT-HR-°F 
D * cylinder, diameter- — Ftr. 


U« * free stream velocity - FT/HR 

p * density - LB/ FT 3 

v « viscosity - LB/FT-HR 

C D v 

Pr - Prandtl number - — dimensionless 

C p = specific heat - BtU/LB-°F 

C ■ empirical constant dotettnined from Table 2-3 
as a function of Reynolds number (Re) ■ U« p D/y 

n » empirical exponent determined from Table 2-3 
as a function of Reynolds number (Re) ■ U*.-p D/y 


OJM 

4-40 

40-4,000 

4,000-40,000 

40,000-400,000 


0>691 
-821 
.615 
.174 
.0230 


!5 



2. 2. 2. 4. 6 Forced Convection Over Spheres 


Empirical relationships, for. average value heat transfer coefficients 
for flow. over-spheres are- given for gasses and- liquid as. outlined in. the 
following steps t 

1) Calculate Reynolds number Re *-U® e_D/y 
where: 

U„ = free stream velocity — FT/HB. ... 

p - density * LB/FT 2 
D ■..sphere diameter - FT 
y ■ viscosity- LB/FT- HR 

2) For a gas where_(T < Re < 25) use step 3 

For a gas where (25 < Re < 100,000) use step 4 
For a liquid where (1 < Re < 2,000 l) use step 5 

3) Gas where 1 < Re < 25 

h = Gp U® p®^ Re + j 

where: 

h ■ convective heat transfer coefficient - BTU/FT 2 -HR-°F“ 
Cp ■ specific heat - BTU/LB^E 
U„ ■ free- stream velocity * ET/HR 
p® «*• free stream-dens1ty_r UB/FT* 

Re » Reynolds number from- Step 1 - dimensionless 

4} Gas where 25 < Re * 100,000- 

H . 8.37 


y 


where: 

h « convective heat transfer coefficient - BTU/FT 2 --HR- 6 F 
k = thermal conducti vttlLr BTU/FT-HR- 0 F 
D * sphere diameter - FT 
U« = free streamvelocity — FT/HR - 

= free stream density - LB/FT! 

u = viscosity..- LB/FT-HR 

5) Liquid where 1 < Re- < 2000 

h * |^b.97 + 0.68 j 

where: 

h = convective heat transfer coefficient - BTU/FT 2 -HR-°F 
k - thermal conductivity - BTU/FT-HR-°F 
Pr = Prandtl nusber = - dimensionless 

Cp = specif 1< heat' - BTU/LB-°F 
u a viscosity - LB/FT-HR 
D ■ sphere diameter - FT 
U* = free stream velocity - FT/HR 
p« » free stream density - LB/FT 3 
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2. 2.2.5 Computational Methods - Radiation Conductors 

Mosi-theanal analyzer: computer-programs linearize the radiation term 
prior to performing the. Heat balance at each time step. This_ope ration 
simply amounts to factoring Jit 4 - Tji 1 ). Into the following components 
(Tt 3 + Ti Tj* + T i. 2 Tj" + J Fj 3 )(Ti - Tj), the. term- (Tj 3 + It Tj 2 + 

Tl z _Tj_* Tj 3 ) is- evaluated by the -computer each time pass using the 
current- values of.Ti and Tj. Thts-quantlty Is then multiplied by the In- 
put- value, of the-radlatlon conductor -thus .reducing the_rad1ation. equation 

to a linear form. The thermal engineer need only be. concerned with the 

Input value of the radiation conductor Which takes the following form: 

G aoeF-j.j At for radiation to a black bedy_ 

G-* for radiation between grey surfaces 

where: 

G » input- value for radiation conductors - BTU/HR-°R 4 
a = Stephan-Boltzmann constant * ,1713xlG' 12 - BTU/HR-FT 2 -°R 4 
e * emittanee - dimensionless 

F|_j ■ geometric- configuration factor from surface 1 to 
surface jr- dimensionless 

At « area of surfaee i - FT 2 

■ grey body radiation factor - dimensionless 

The emittanee, «,1s a measure of how-well a body can. radiate energy as 
comparod with a black body. Emittanee is the ratio of the total-emissive 
pew?f of a real' surface at temperature T te the total emissive power of a- 
blaek surface at the same temperature. The emlttanees of various surfaces 
are- a function of the material, surface condition, and temperature of a 
body. The- surf ace of a- body, and therefore the- emittanee, may bealtered- 
by polishing* roughing*, painting, etc; The values of e for mest-eonmon 
materials and surface conditions- have been- measured at Various tempera- 
tures and are presented In tables or graphs in many reference manuals. 

It is left te the engineer to determine the value of enrfttance to Be used 
and-whether the variation of t with temperature- Is significant oyer the 
temperature range expected for the surface. 
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The geometric shape (configuration) factor from surface 1 to sur- 
face j, Fi-j, Is the fraction of total radiated energy from- surface 1 
which Is directly Incident on surface i, assuming surface 1 to be emitting 

energy. diffusely . Fj_-j would be the fraction of' total radiant. energy. 

from surface j- which Is- Intercepted: by surface 1. The configuration 
factors for finite regions. Of jdlffi use -areas are related by: 

A i F 1.J * A j F j-1 

The configuration facto r,_F^_^, Is a function of the geometry of the system 
only. Several . iter programs have teen developed_to compute the shape 
factors between surfaces with complex geometries! however, form factors 
between some surfaces with simple geometries can be hand, computed. 

Figures 2.22 through 2.38 present configuration factors for various 
simple geometries. The use of these Figures and-conflguratltm factor 
algebra will allow the engineer to determine form factors for many simple 
radiation problems'. 

The following examples of configuration factor algebra shouldube 
helpful; 



A l-2 F 12-34 " A l F i-3 + A l F l-4 * * ^2-4 

A|F^_ ; * A^F-j^ (symmetrically positioned} 



Figure 2-22. Radiant- Interchange Configuration Factors Paint Sources 
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This configuration is a plane point source «Mj and a plane rectangle A% 
parallel to the pl*ne of dAj. The normal to dAj passes through one-cor- 
ner of As. The curves for this configuration are given above where Fjg 
Is plotted as a function of * and y,~w1th x « a/e and y * b/e. 
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This configuration Is a spheric*) point source <Mj and a pi an* rectangle 
As 5 the point srurce Is located- at one corner of a rectangle that hat- one 
common side with As. The planes of the Wo rectangles Intersect at an an- 
gle a. The configuration factor Fjs Is plotted above as a function of 
x and y, where * <• 5/a and y • a/o. 




from NASA Report No.- TN-4B36,- by D. C. Hamilton and W. ft. Morgan 
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Figure 2-22 (Continued) 



This* configuration Is a pUr> mint source dAj and any Infinite- pi ane-Ag, 
with the planes of dAi and Ag Intersecting at an angle e. The configuration 
factor values may be calculated from tbs above equation. 



? t2 — —(cos 6 — cot-G>> 


This configuration is a plane point source dAi and any surface Ao generated by 
an Infinitely long line, moving parallel te itself and to the plane of dAj. 
The configuration- factor values may be computed from the above equation. 
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Figure 2-22 (Continued) 
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This configuration Is a plane point source dAj and a plane circular disk 
42 . The plane of dAi Is parallel to toe plane of An the point source i: 
located at a distance n from the normal to the center of 42. The conflgu* 
ration factor Pjz Is plotted above for various values of £ and D, where 
s • rz/d and D ■» d.'?2- 
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TiHs configuration Is a plane point source dAi and a plane disk *21 the 
planes of tttj and 42 trrtersect at an angle of do°> The- centers of 42 and 
<24 lie In a plane that Is perpendicular to the two planes. The equation 
above Is given 1ft terms of R and P for R »< r/b and P « a/b . 



Figure 2-22 (Continued) 
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This configuration Is a plane point source- d<4 and a right circular cylinder 
As of length- 1. The normal to dAi passes through the center of one end of the 
cylinder and Is perpendicular to the axis of the cylinder. The configuration 
factor Is plotted above as a function of D and L where D - d/r and L - z/r. 



I m ’ m l,,M V 1 » » ■•* • 
lMSa^>^ilMkk!'IVIk>TPW%;l 

■n r ^!tiik?vn^iaG*iiil 

isr^a::Vit»?wk^ 




Ml 



1111 




ew^tasl 

/IIIMK«kXk.'lkM 






3 4 9 7 10 (9 20 3040 
D 


This conffguratlon consists Of tiro concentric cylinders of radius r and d 
and length l with a point source dAj on the Inside of the la ge cylinder at 
one end. The configuration factor Pji Is from the point source <Uj on Ai to 
AliAl does not Include the ends of tne annulus. The curves for this configu- 
ration are given above In term? nf d and i, with D ■ d/v and £ » l/r. 










Figure 2-22 (Concluded) 



This configuration consists of an Infinitely long cylinder Aj and an 
Infinite plane Ab, mutually parallel . M - m/r and N » n/r. 




Figure 2-23, Radtafit Interchange Configuration. Pastors 
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This configuration Is a line source dAi and a plane rectangle A% parallel 
to the plane of dAj t with dAj opposite one edge of As. The plot of Fjs Is 
a function of * and y. where x * S/d and y • a/a. 
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This .configuration is a-Hne source d4i and a right circular cylinder /ig, 
both of ' .length- Z; the no mat -through each-end of the source- passes through-, 
and'lt 'normal to, the center line of the cylinder at- the ends. The configu- 
ration factor values ar* 7 plotto<f above- In terms of the parameter s and s, 
where D » d/p and~£ » I/p. 


f Adapted fi*em NAGA Repert-N<^'N-4B36* by C. HanrtRon and W, R. Morgan 





Figure 2-23 (Concluded} 
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The cr<*y-oody shape factor (Script F )^j» Is the preduct-of the 
geometric shape facte r — F-j„j and_a factor which, el lows for the departure of* 
the surface from black, body conditions* For= radiation enclosures, the 
^ j factors are generally evaluated with a computer program. The input 
for the program being the A i F j-_j values, for every surface of toe enclosure 
to every other surface and. the emittance and area -foreach sur.face, __S1mpl i ffed . 
e quations for.j^ ^ exist'for two-componentrgray enclosures. 

Parallel flat, plates;- F 1 ^ « F^ = 1 


Concentric cylinders of Infinite height or concentri c spheres : 
F-j ,2 “ i» p£-l ^ ® 



For "non-enelosed M surfaces an effective emittance,. th* 

surfaces may be used to compute the gray-body form factor with toe following 
equation: 


•^1-j “ 6 eff_l F 1-j 

The effective .emittance is a. functlon-of the-emtttanees of the- Re- 
surfaces and toe- configuration factors (1) be fcween-them. - Tfie error Induced 
with use Is the-result of -neglecting secondary reflections from 

surfaces other than the two for which the- effective- Omittance was determined. 
By-reduci ng-Hottel* s method for three flat plate surfaces with ami ssfvi ties 
of i-\ t eg:, and 1« the following equation ean be constructed: 

« e l e 2~ 

eFT 1 * F j o 1 ^ “ e l^ ^ ^ 


Tabic 2-4 presents approximations for' various cj>, and 
f 1-2 * p 2-t v * JueS ‘ Intermediate values can be approximated by- interpolation. 
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Table 2-4 (Continued) 


flat plates inside a slack enclosure 
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.263 

.225- 

-.!«1 

.130 

.070 

.017 

.008 

0 

.30 

*100 

-.286 

.270 

.232 

.232 

.208 

• IP* 

.168 

.109 

.080 

.032 

.007— 

F 

.30 

.200 

.192 

.183 

.173 

.161 

.167 

.136 

• 109 

.082 

• 067 

.026 

.005 

C 

.10 

.100 

.097 

.093 

.089 

.086 

.078 

•070 

.060 

• 047 

.028 

.016 

.003 


.06 

.060 

.069 

.067 

.065 

.043 

.060 

.037 

.032 

.026 

.016 

.009 

.002 


.01 

.010 

.010 

.010 

.009 

• 0o9 

.008 

.038 

.007 

• 005 

.003 

.002 

>000 


FLAT PLATES INSIDE A SLACK ENCLOSURE 
EFFECTIVE EMITfAftCE FOR F12«F2l • .70 


<»HTT'ANCE OF SURFACE 2 




1.00 

.90 

.80 

.76 

.60 

.30 

.60 

.30 

,20 

.10 

.03 

.01 

t 

M 

1.00 

1.000 

• 900 

•son 

.700 

.600 

• 300 

.600 

• 300 

.200 

•lOO 

•030 

•010 

I 

.90 

.900 

.816 

<730 

.6*6 

.556 

.666 

.176 

.286 

.191 

«096 

• 068 

• 010 

T 

t 

.80 

.800 

.730 

.638 

.385 

.306 

.630 

.369 

.266 

.180 

.092 

• 6«* 

.009 

A 

N 

f 

.70 

.700 

• 666 

.$83 

.323 

.639 

.391 

.320 

.266 

.168 

;o86- 

.046 

.008 

I 

.60 

.600 

.636 

.508 

• 639 

.664 

.369 

.268* 

.226 

.133 

.060 

.061 

• 009 

0 

.30 

.500 

1666- 

.630 

.391 

.369 

;303 

.253 

.199 

.m 

.073 

♦037 

.008 

P 

.60 

• 600 

.376 

.349 

*420 

iitft 

• 233 

.216 

.170 

.120 

• 066 

.033 

<007 

$ 

u 

o 

.30 

.300 

.286' 

• 266 

.268 

.228 

.199 

.no 

.117 

• 099 

• 036 

.028 

.006 

R 

F 

.20 

.200 

.191 

.189 

.161' 

'.153 

• 139 

.129 

.099 

.072 

.060 

,<ni 

.006 

• 

c 

• 10 

>100 

.096 

.092 

.081 

.080 

•dls 

.066 

.036 

• 060 

.023 

.012 

•on! 

E 

1 

*03 

.050 

.068 

.nV6 

*066 

.061 

.037 

.033 

.028 

• 021 

*0i2 

.007 

OOl 

• 01 

.010 

.010 

.009 

.009" 

”1008 

•008 

.out 

.006 

.006 

.003 

•0O1 

•non 




Table 2-4 (Continued) 


flat plates inside a black enclosure 

EFFECTIVE tNITfANCe FOA 312**21 m , 6o 


CNITTANCE OF SURFACE— 2 




1.00 

.90 - 

.80 

*20 - 

-.60 

• 60 

.40 

- .30- 

.20- 

.10 

.06 

• 01 

E 

M 

1.00 

1.000- 

*900 

.800 

• TOO 

.600 

• 900 

*400 

.300 

.200 

.100- 

.060 

•010 

1 

f 

.*0 

.900 

.816 

.729 

• 642 

•963 

.464- 

.379 — *292- 

- .189 — .096-. 

-.640 

• Ole 

T 

A 

.80 

.900 

.729 

.636 

*681. 

•604 

.426- 

*346 

.262 — .177 

•090 

.046 

.009 

N 

C 

.TO 

• TOO 

.642 . 

.581 

.518- 

- .463- 

*386 — .614 

• 240 

. 164 

.084 

.042 

.009 

E 

.60 

.600 

.553 

.904 

.463 

.398 

<341 ,280 

.216 

• 149 

.077 

•039 

.008 

0- 

F 

.SO 

.600 

.444- 

.-426 

.396 

.341— .294 

.244 

.190 

.132 

.068 

.035 

.007 

S 

.AO 

.AOO 

• 373 

.349 

<314 

.280 

.244 

.204 

.160 

.112 

•059- 

.030 

.006 

U 

R 

.30 

.300 

.282 

.26- 

.240 

.216- 

.190 

. 160 

.127 

• 0*0 

.048 

.026 

• 005 

F 

A 

.20 

.200 

.199 

.177 

.164 

.149 

.132 

.112 

.090 

.066 

•036 

.018 

.004 

C 

E 

.10 

• 100 

.096 

•09n 

.084 

.077 

.068 

.069 

.048 

.036 

♦ 019 

.010 

.062 

1 

.06 - 

--.OSD 

.048 

.049 

.042 

.0)9 

.036 

.060 

.026 

.018 

•olo 

.006 

• 001 


.01 

.010 

♦ 010 

.009 

• 009 

.008 

.007 

.006 

• 006 

• 004 

•002 

• 001 

>000 


FLAT PLATES INSIDE A BLACK ENCLOSURE 

effective enittancE FOR Pil'fn m .SO 


WIT?a*|C£ OF SURFACE "2 




1.00 

.90 .80 

*70 

.60 

.60 

• 40 

*10 

.20 

*10 

.06 

.01- 

e 

M 

t.ao 

i.ooO 

•900 .son 

• TOO 

•600~ 

O 

o 

21 

«. 

.400 

• 300 

.260 

.100 

•040 

.010 

r 

T 

.90 

.900 

.814 .727 

.640 

•661 

.462 

.371 

• 2io 

.187 

.094- 

.047 

>009 

T 

A 

• 80 

.800 

.727 .413 

.577 

.500 

• 421 

• 340 

.266 

.174 

• 066 

*044- 

•609 

N 

C 

.70 

.Toft 

.*40 .977- 

• 613 

.447 

.378 

• aOS 

,»66 

.169 

.0(1 

.041 

.006 

E 

.60 

.600 

•561 .500 

.647 

4361 

.323 

.213 

.209- 

• 149 

.073 

.037 

• 607 

0 

F 

.50 

.600 

• 462- .421 

.3*6 

.333 

.2(6 

.23* 

.1(2 

.125 

.666 

.6)3 

.007 

S 

.40 

.400 

•371 *340 

• 108 

.273 

• 236 

.196 

.162 

• 109 

.066 

.028 

4006 

u 

R 

.30 

.300 

• 260 •»! 

.221 

.26* 

• 162 

.192 

.119 

.083 

• 064 

.022 

•60S 

F 

A 

•20 

.200 

.141 .174 

.199 

.143 

*126 

.106 

• 063 

• 066 

.031 

.018 

.00) 

C 

c 

.10 

.100 

.094 , 04r 

• Oftl 

•073 

.066 

.095 

• 044 

.061 

• 017 

• 009 

• 002 

1 

.06 

.060 

.047 *044 

.041 

.037 

.031 

.029 

•022 

.014 

• 009 

•006 

.001 


.01 

.010 

.009 .009 

*008 

.007 — 1007 

• 006 

.009 

.003 

.002 

.001 

.000 
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Table 2-4 (Continued) 


FLAT PLATES INSlOt * BLACk ENCLOSURE 
EFFECTIVE EHtTTANCE FOR F1Z*F21 • .AO 


EHITTANCE-JOF surface 2 



1 .nn 

. Qft 

ft IV. 

.70 

*60 

.30 

.40 

.30 

.20 

. 10 

*03.- 

.*01 


*• 00 

• "O 

• Bv 




1.00 — 

l.ooo 

.$00 

.806 

.TOO 

.600 

.300 

.400- 

.300 

.200 - 

.100 

*030 

•010 

.90 

•«00 

.813 

.726- 

*638 

.349- 

.459 

.369 

.278 

.186 

• 093 

.047 

.009 

.80 

.800 

.726- 

*630 

«S2A 

*496 

• 4U- 

.336 

.254 

.171 

• 086 

.043 

.009 

.TO- 

. .TOO 

.638 

.526 

.308 

.44 1. 

.372 

. 302— 

.229 

.153- 

.078 

- .040 

.008 

♦60- 

.600 

.34$ 

.496 

-.441 

.383 

*326 

.263 

.203 

.138 

.070 

-.033 

.007 

.so 

.500 

• 43V 

.417 

.372 

.326 

.218 

.227 

.174 

.119 

.061 

.031 

.006 

.40 

.400 

.369 

.336 

.302 

.263 

•22J 

.187 

.144 

.099 

.051 

.026 

.005 

.SO 

.300 

.278 

.234 

.229 

.203 

.17# 

.144 

.112 

.077 

.040 

.020 

• 004 

.20 

.200 

.186 

.171 

.133 

.138 

.119 

.099 

.077 

.034 

.028 

.014 

.003- - 

.10 

.100 

.093 

• 0<6 

•078 

.070 

.061 

.051 

.040 

.028 

.013 

.008 

.002 

.03 

.030 

.047 

.043 

.040 

.033 

.031 

.026 

.020 

.014 

•ooi 

.004 

.001 

.01 

.010 

.009 

.009 

.008 

.007 

.006 

.003 

.004 

.003 

.002 

.001 

.000 


FLAT I'LATES INS tot A BLACK ENCLOSURE 
EFFECTIVE ENITTANCE F6h F12*F21 • .30 


EMI If ANCE OF SURFACE 2 



i 

.00 

.90 - 

.*80 

-.70 

.60 

*50 

.40 

.30 

.20 

•10 

•os 

.01 

E 

U 0 O 1 

*000 

.900 

.800 

*700 

*600 

.300 

.406 

.300 

• 200 

• ids 

1030 

.010 

H 

1 

.90 

.900 

.812 

.724 

1636 

.347 

*437 

.367 

.276 

.164 

.092 

.066 

•009 

T 

T 

.80 

.800 

.724 

.648 

.370 

.46* 

.412 

.132 

.231 


•042 

.009 

A 

N 

.70 

.700 

*636 

.370 

.504 

.436 

.366- 

.296 

.224 

.131 

.076 


C 

E 

.60 

• 6 o 0 

.447 

.4*2 

.4f6 

.376 

.318 — .259 

.197 

.133 .('67 

• 036 

.007 

0 

.50 

.300 

.437 

.*» 

.766 

.319 

.270- 

# 22 o 

• 166 

.114 

.938 

.029 

.006 

F 

.40 

.400 

.367 

.332 

.296 

.239 

.220 

.179 

.m 

.093 .046 

.024 

.003 

& 

U 

.30 

.300 

.276 

.231 

.2*4 

*197 

.168 

.117 

*106 

*072 

.037 

.019 

,004 

ft 

F 

.20 

.200 

.184 

.168' 

.131 

.133 

.114 

10*3 

*072 

.030 

.026 

.013 

.003 

A 

C 

.10 

■1100 

*092 

.0C3 

.076 

.067 

*038 

.640 

*037 

.026 .013 

.007 

.001 

E 

.05 

*030 

.046 

.042 

•038 

.03* 

.029 

.024 

.619 

•oi i 

♦ 007 

.003 

•oof 

1 

.01 

.010 

.009 

.004 

• 008 

.587 

.006 

.003 

• 004 

.003 

.not 

.001 

• nnn 
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Table- 2-4 (Continued) 


flat plates inside a black enclosure 
EFFECTIVE ENltTANCfc FOR F12«F21» .20 


EMJTTANCE OF SURFACE 2 




1.00 

.90 

«00 

.70... 

.60 

.50- 

-“•40 

.30 

.20 

*10- 

♦08- 

-.01 

£ 

M 

uoo 

1.000- 

.900 ... 

-.900 

.700 

•600 

.500 

.600 

.300 

.200- 

.100 

• 050 

*010 

L 

J 

*90 

.900 

.012- 

.722- 

.636 

• 544 

.655 

-.366 

.276 

-.183- 

-.092 

-•066 

• 009 

T 

A 

.to 

.000 

.723 

• 664 

• S67 

.688 

.608- 

.328 

.267 

• 165 

.083 

.062 

>008 

H 

c 

*tg 

• TOO 

636 

.567 

*699.. 

.630 

.361 

.290 

.219 

.167 

.076 

.037 

.007 

E 

•to 

.600 

.466 

*608 

.630 

.372 

.312 

.252 

.191 

.126 

.065 

.032 

.007 

0 

r 

.so 

.100 

.655 

.600 

.361 

.312 

.263 

.213 

.161 

.109 

.055 

.028 

.006 


A.*0 

.600 

.366 

.320 

.290 

.252 

.213 

.172 

.131 

.088 

.063 

.023 

.005 

u 

.30 

.300 

.276 

• 24? 

.219 

.191 

.161 

.131 

.100 

• 0&ft 

.036 

.017 

.003 

p 

1 

.20 

.200 

.103 

.165 

.167 

.128 

.109 

•008 

.068 

.066 

.023 

.012 

.002 

c 

£ 

.10 

.100 

.092 

.003 

.076 

.065 

.055 

.065 

.036 

.023 

.012 

*006 

.001 

1 

.OS 

.oso 

.066 

.062 

.037 

.032 

.028 

.023 

.017 

.012 

• 006 

.003 

.001 


.01 

.010 

.009 

• 008 

• 007 

.007 

.006 

• 005 

• 003 

.002 

.001 

.001 

• OOO 






FLA7 

PLATES 

INSIDE 

: A SLACK ENCLOSURE 








EFFECTIVE EM1TTANCE 

: FOR F 12» F 21 

. • • 

10 









EMltTANCE-OF SURFACE 

2 







1.00 

• 90 

.89 

.70 

.60 

.90 

.60 

.30 

• 20 

.10 

.03 

•01 

E 

M- 

1.00 

1.000 

.9 no 

• 80fl 

.700 

.600 

.900 

.600 

.300 

• 200 

• IOO 

.090 

>Alo 

I 

T 

.90 

.900 

•an 


.632 

.562 

.692 

• 162 

.272 

.im 

.091 

.043 

.009 

T 

A 

.00 

.800 

.721 

.061 

.963 

.666 

.606 

.326 

.263 

.163 

.081- 

• 041 

.008 

N 

r 

.TO 

.700 

.6 32 

.961 

.696 

.625 

.395 

.209 

.219 

.163 

.072 

.036 

.007 

i 

.60 

.600 

.162 

.666 

.625 

.366 

• 906 

• 266 

.189 

.124 

.062 

.031 

• 006 

0 

p 

.30 

.900 

•6$ 2 

.606 

.339 

• 306 

.296 

.206 

.1)3 

<l66 

*&32 

.026 

• 005 

4 

• 60 

.600 

.312 

.326 

.265 

.266 

.206 

• 166 

.124 

.066 

• 062 

.021 

• 004 

u 

R 

.30 

.300 

.272 

.26) 

.213 

.163 

.193 

• 129 

• 099 

• 066 

•012 

.016 

•oot 

F 

A 

.20 

.200 

.111 

.109 

.14) 

.124 

.106 

.086 

.086 

.063 

• 022 

•Oil 

.002 

C 

£ 

.16 

.100 

• 09 i 

.081 

.072 

•062 

.032 

.062 

«C 12 

.022 

.011 

.005 

.001 

1 

.03 

• 090 

• 069 

• 06] 

• 036 

.031 

.026 

• 021 

.016 

.Oil 

• 003 

.00) 

.001 


.01 

.010 

.009 

.006 

.007 

• 006 

.003 

.006 

.001 

.002 

.001 

• 00) 

.000 
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Table 2*4 (Concluded) 


FLAT PLATES IflSlDt A SLACK- ENCLOSURE 
EFFECTIVE EMITTANCE FOR F12*F21 • .06 


ENIITAKCE OF SURFACE 2 



1.00 

.90 

.00 

.TO- 

.60-- 

E 1.00 
M 

1.000 

.900 

.600 

•TOn 

.600 

I .90 

y 

.900 

.010— .321 

.621 — 

.641 

T .SO 
A 

iSOO 

.721 

.641 

• 962 

.612 

K .TO 

.TOO 

.631 

.662- 

.692 — .623 

C 

e .*o 

.600 

.541 

.402 

.423 

.363 

o .so 

c 

.600 

.491 

.402 

.393 

.303 

.40 

.400 

.361 

.322 

.203 

.243 

S 

U .30 
A 

.300 

.271 

.242 

.212 

.163 

F .20 
A 

.200 

.101 

.161 

.142 

.122 

C .10 
E 

.100 

.090 

.001 

• 071 

.061 

.03 

9 

.060 

.049 

• 040 

.036— .031 

.01 

.010 

.009 

• 000 

.007 

.006 


.30- 

.40 - 

-.10 

• 20 

.10 

.09 

.01 

.300 

.400 

.300 

.200 

•mo 

.090- 

.010 

.491- 

-.361 

.271 

.161- 

.090 

.069 

.009 

.402 

-.322- 

.242 

.161 

.001 

.040 

.000 

.331- 

.201- 

.212 

.162 

.071 

.036 

•007 

.303 

.243 

.133 

.122 

• 001 

• 031 

.006 

.233 

.203 

.133 

.102 

•OM 

.026 

.009 

.203 

.163 

.123 

.002- 

.041 

.021- 

.004 

.133 

*123 

.092— 

.062 

•031 

.016 

•003 

• 102 

.002 

•062 

•041 

.021 

.010 

.002 

.031 

.041 

.011 

.021 

.010 

.003 

.001 

.026 

.021 

.016 

.010 

*003 

•003 

.001 

.009 

•004 

.003 

.002 

•001 

• 001 

.ono 


Ft AT PLATES JNSIOE A SLACK ENCLOSURE 
EFFECTIVE EHITTaNCE FOR FI2*F21 * .01 


CMlTTAneE OF- SURFACE 2 


e 

1.00 

1.000 

.900 

*600 

• Too 

•too 

.906 

.400 

1300 

.200 

.100 

.030 

.010 

t 

f 

.90 

.900 

.110 

.720 

.030 

.940 

.430 

.300 

.270 

.100 

.690 

.043 

.609 

t 

A 

.30 

• 

Q< 

O 

o 

.120 

.640 

.300 

.436 

.600 

.320 

.240 

.160 

.000 

•040 

.00# 

H 

.70 

.TOO 

.630 

.36ft 

.490 

.42* 

.331 

.261- 

.210 

.140 

.070 

•039 

.06? 

e 

.60 

.600 

.340 

•4.3 

• 421 

.361 

.301 

.241 

.m 

.126 

.060 

.0)0 

.004 

0 

r 

.SO 

.600 

.430 

.400 

.391 

.101 

.231 

.201 

.131 

.160 

• 090 

.029 

1009 

« 

.40 

.400 

.160 

• Ml) 

.331 

i»l 

.201 

.161 

.121 

2000 

•040 

.020 

.004 

u 

U 

■ 30 

*300 

.170 

.240 

.210 

• Ml 

.1*1 

.131 

.090 

.060 

.6)0 

•019 

.009 

r 

A 

• 20 

.200 

.lio 

• too 

*140 

.120 

1100 

*630 

.060 

.040 

• 020 

•oio 

.662 

n 

c 

• 10 

.100 

.090 

.0*0 

.070 

.040 

•030 

.044 

.010 

.020 

.010 

.003 

•001 


• 03 

.0*0 

.043 

.046 

.039 

.030 

•023 

•o3o 

.019 

• 010 

.063 

.003 

.001 

a 

.01 

.CIO 

.069 

•oot 

.007 

.006 

.003 

.004 

.003 

• 002 

.001 

.001 

.000 
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2. 2. 2. 6 Computational Methods - Mass Flow Conductors 


The use of a mass_flow-conductor In a thermal network Is a convenient 
method of accounting for the transfer of energy, from one point to another 

due to the actual movement fflow) of. a fluid from one point to another. 

Mass- -flow conductors are computed from the equation: 


wherer 

* = the mass flow rate of. the fluid.*- LB/HR 
C p « fluid specific heat — BTU/LB-°F 

The mass flow rate {*) Is related to the fluid velocity by the expression: 


W 3 p A U 

where: 

A - the cross sectional area through which the 
fluid flows - FT* 

U « the fluid velocity - Ff/HR 

p 3 the fl4ii d densl - LB/FT? 

The thermophysical properties Cp and p should be evaluated at the bulk 

temperature of_the fluid. 

The mass flow- conductor which Is In addition to any other mode of- 
heat transfer sim&ly accounts for the Internal energy term of a mass 
moving from one location to another. 



2.2.3 Energy Sources or Sinks 
2,2. 3.1 Concepts 


Energy sources or sinks, Q, ar'e modeling elements which allow. the 

Impression of positive or negative- heating rates on the nodes of a thermal- 

network Independent of. conductor paths to the node. ... 

2.2. 3. 2 -Types of Heat Sources or Sinks 

Common engineering applications of heat sources In thermal models - 

are: 

• Solar and Planetary_Heat1ng 

• -Aerodynamic Heating 

• Avionic Col dpi ate Heat Loads 

• Change of State Latent Energy 

• Thermal Control Heaters 

Common application for heat sinks are: 

• Change of State Latent Energy 

• Radiator Heat Rejection 

• Aerodynamic Cooling 

2.2. 3. J Computational Considerations - Sources or Sinks 

Heating rates_may be impressed on diffusion (finite capacitance} or 
arithmetic (zero capacitance) nodes. Most thermal analyzers provide a 
separate entry block for entering .heating or cooling rates. For example, 
the SINDA computer- program uses the SfiURCt DATA BLOCK for such entries. 

In the usual case, heating rates are not considered when computing the 
time steps for transient analysis, and large heatingirates on low capaci- 
tance nodes mgy. create instability In the network sotutlon. Also the 
impression of large heat sources on arithmetic nodes with radiation (non- 
linear) conductors- attached often Causes large erroneous temperature 
oscillations In the arithmetic and adjoining nodes. Both of these 
difficulties can be avoided with the use of the program control constants In- 
corporated in most thermal network analyzers. These control constants are 
the time step multiplication factor and the maximum temperature orange allowed. 
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2.3 Network Solution 

An area which Is of utmost concern to the. engineer Is an-explanation 
of the system of equations which are solved in a typical thermal problem. 

At this point it must be completely understood that- the physical system is 
reduced to a lumped_(nodalized) system and that: the_choi ce of the nodal 1- 
Zattomhas a far greater bearing on the problem than -the parti eutar 
numerical -solution-technique and asrsumpttons—that are usually_employed in 
the problem solution*- Each numerical solution me thocLis~ bounded by 
accuracy, constraints-dictated-by the_pacticulac_set of assumptions: assa- 
ciated with the particular problem undergoing solution. 

The basic transient heat transfer equation applied to linear conduction 
problems isr 

* 9* T ,__2? . ^L+ *L + ,,, 

9t 3X 2 ay 2 az 2 ( » 

The solution techniques most commonly applied stem from a reduction of 
Equation (1). It can be easily seen that any nonlinear terms such as 
radiation must be linearized by some method if the above equation Is to be 
used. 

Eor example the radiation term GOi 4 - Tj 4 ) can be linearized by 

G Y (T 1 - Tj) - G(Tl 4 - Tj 4 L (E) 

where-y is tbe^ linearization factor. The- calculation- of y and the lineari- 
zation of the G^s-is performed automatically by thermal network analyzersL 
and will be explained in more detail in the- forward differencing. section. — 

The-pracess of nodaTization reduces the- volumetric dimensions and- 
properties so that equation (1) »‘an be cast in the following femn- 
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Each tens of the summation, ®-jj(Tj-Tj)> represents the- heat rate^jQ. 
following Into node - 1 from node - j, -Integrating equation (3) yields: 



To actually, pejcformitbe^i-ntegratlon Jndi cated in equation- (4-). on a computer 
It would' be- most' convenient -If- the.. integrand war not a function of time: 


tnew __ 

-T.| ) | = = constant — (5) 

*old- - — — — — 

. If this was the case, then the second term in equation (4) would reduce tor 
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old 


new 
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W T i 


) dt 


Unew * told) 
— 


'1 




at. 


( 6 ) 


and equation (4} could be east in the simple^ form of: 

Qi 

T^fc* at) « T^(t)> ^ fit- (Z)- 

Equation (Z). represents the basie finite- differencing formulation of the 
transtentheat transfer equation. AIT finite differencing methods-, be they 
forward, backward, central, explicit or- tmpTI ett^_perfom the calculation 
1ndtcated-4n equation (71. The only-difference among- the various Emula- 
tions Hes-ln the- assumption used- to evaluate the- Qi term-per equation 
The significance of this assumption as it relates to the various finite 
dITf e rcr.ci ng_fo tmulatlons Will- be discussed in following seetiohs. 


2.3.1 Stead«_State 


A thermal system has_reached the steady state when the net heat flow- 
rate to each diffusion jnd arithmetic node fs zero. From equation- (3), 
this yields: 


\ * °- 


( 8 ) 


Rearranging the terms In equat1qn(8), yields 


V 


Ml 




(9) 


The subscripts In equation .(9) obscure the- elegance of the steady state 
solution, which can also be expressed as: 


I ■ [63 I 00) 

Equation (10) clearly reveals that finding the steady state temperature 
vector simply- amounts to finding the eigenvector of [G] for an eigenvalue- 
of 1, The only problem is that, for a- network with more than a few (say 
20) nodes ,.the-ca1 dilation of- the^ei gen vector using an explic-ltnalgorlthm, 

such as. Gaussian Elimination, would take a prohibitive amount-of time 

this t. in addition, to the fact that-most of-the elements of— [G] are zero,, 
suggests that - an iterative- technique such as expressed in e quation (U) 
would result in a. huge reduction- In. requl red computer time. 

f t+1 - m T- 1 " (H) 

When It is also realized that-equatlen- (18) is 6nly-val44-when the G-mafcrlk 
Is absolutely constant with- respect to temperature* it becomes clear* that, 
equation (11) is: the only solution technique that is universally usabie-j 
Using tills technique* it Is simple-matter to update the G matrix (evaluate 
temperature varying G r s and linearize radiation 6's) based oh the new~gi‘ess 
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ttmiperature— vector T,^ 4 ^ before-proceeding to the next iteration.- In 
addition,, equation (11)_ does_not require- that the-T 1 ^ and T^ectors be the 
same- length, so that the-jJ vector, can he extended to Includes vector of 
temperatures which represent the boundary conditions placed on the network. 

Two criteria for terminating_the- iterative .process are generally 
provided. The ftrst-is simply a-flxed number, N, which represents the — 
maximum-number of iterations to be performed- The second Is a convargence- 
ar, as It is-moce_commonly called,, a, relaxation criterion, 6-, which is 
defined asr : 

max Tj ^ - • - <_ & Cl2) 

j 

In other words, calculations cease when no "significant" Improvement from 
one Iteration to the next Is noted (where "significant" is defined In 
terms of s). 

Strong oscillation In the temperature vector Is often noted when 
analyzing a network which Is dominated by radiation effects. To hasten 
(and in some cases, to enable) the convergence of the solution* a damping 
factor, c * Is often applied to equation (11) to yield the more versatile 
solution technique expressed In equation (13>. 

t 1+1 - C [G] i* + (l-sl-I 1 

to < ctl) (U) 

To summarize, the fully general finite differencing technique for 
steady state analysis is- an iterative algorithm requiring three parameters, 
N, S', and $, for the control of convergence. It might be noted here that 
arithmetic nodes always receive a steady state solution technique even 
when they appear In a -network which Is undergoing a transient analysis* 


D 
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O 2.3.2 Transient Analysis 

2.3*2. 1 Forward Differencing 

Forward differencing derives Its name from the fact that all tempera- 
tures are extrapolated forward In time, t,_for the purpose of evaluating the 
expression in equation (5). This point is illustrated in Figure- 2-26. 



Figure- 2-25. Forward Extrapolation Used in Forward Differencing 

t "i) Combining equations (5) and (7) under the above assumption yields the basic 

forward differencing equation: 

T,(t*4t) - T,ft) +gi £ G^ft) - V‘)! <14) 

1 i 

As might be expected, the temperatures at tlme-t are used- to evaluate any 
temperature varying. or. radiation conductor^ That Is, for example* the 
linearization factor, y» is computed as: 

Y» (T^t) + Tf(t>) (Tj(t)- + T^(t)) (15) 

The basic forward extrapolation assumption! has two consequences: 

(a) the resulting solution equation is explicit, and (b)- the solution 
can be unstable. When-it is stated that an- equation is explicit* It means 
that all of the unknowns- are on the' left hand side-and all of the known 
quantities ate on the right hand side-. Since all or the temperatures 
' \ at timet are known, equation (14) defines the temperatures at- time 
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t + At* explicitly. Instability results, from the fact that equation (14) 
is, so to speak* '‘openrended'*. By. choo$1ng_a sufficiently large At, the 
new-tenperatures can be made, unreasonable. . It has been shewn (reference 
3)*. however, that if At is restricted to be less .than, the stability 
factor, t, defined In equation. (15),. then the Solution normally will be 
stable. 

— i_«-C$9UN (16) 

i 1J 

This t, commonly called lithe. CSGMIN*!^ is the smallest time constant in_the 

thermal network. The "CSGMlN" does not include boundary effects. Sharp 
gradients coupled with high magnitude heating rates could possibly cause 
Instability. 

It - is clear, then, that forward differencing represents a two-edged 
sword. On the one hand, the solution equation is explicit, which means 
thatrthe volume of-eomputatlons for each At step will be a minimum. But * 

on the other hand, the size of the time step is limited to the smallest 
time constant- In the entire network:. However, additional advantages 
aecrue from both of these conditions. First* the user need not specify 
any convergence criteria such as were required by lire steady state tech- 
nique, and second, .for normal’ cases the user^need not specify the time step. 

At, since this cart' be conveniently computed by the computer* using a simple 
algorithm such as given In equation (17K 

At * 95* T (17) 

Another moderating point to be considered is that the maximum Atr 
used for a transient analysis Is often-not restricted by the network, 

But rather by the frequency spectrum of-the boundary conditions, m 
gross terms, one cannot use a At of one hour to analyze a network sub- 
jected to a boundary condition which varies at the- rate of 50 cycles per 
hour. In finer* terms. Shannon's Sampling Theorem (rcference-4) dictates 
that: ) 


» s > 2« c 


(18) 



f 


or: 


max At < 


1 

sr 

C- 
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where w is theJilghest: frequency component in the driving signal. Of 
course. ^the limit expressed In equation (19) is theoretical In Jiiat it 
assumes that the- system J s_jl perfect low-pass filter* Thls-ls true, in _ 

practice, only at- nodal potots rather- far away from the boundary. Hence, 
analyses' of. the temperature-response near rapidly varying boundaries: 
will require a At much smaller than the maximum expressed in equation (15). 


It Is no accident, then, that forward differencing formulations enjoy 
the widest and most frequent use In engineering practice. —Other formula- 
tions have practical value only when the choice of At is not restricted 


by boundary conditions 


2. 3. 2. 2 Backward Differencing 

Backward differencing takes Its name from the fact that all tempera- 
tures at time tnew are extrapol ated backward In time In order to evaluate 
the 0 terns per equation-1-15* This point Is illustrated in Figure 2-26. 




Figure 2-26. Backward Extrapolation Used in Backward Differencing 
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Using this assumption- to combine-equations (5) and (7), and stating the 
result-in matrix form yields: 

T-(t+At) » T (ti>-At [G/C3 T (t+At) (20) 

The solution to equation (2) obviously requires that au matrix, inversion 
be performed* but- ft also requires, strictly*, that the G/X. matrl x be. 
temperature Invariant. (Otherwise ,_the basic assumption would require 
that [G/C] be__evaluated for £ « Tj(t±At),_and* of course, .these tempera- 
tures are. not .yet: known.) In. practice*, then, equation. (20)- Is reworked 
into a form which lends i tself to iterative- calculations 


(t+A) 


IT T 1 <*> ♦ 5 

~q 

IT + j S 1 j 


( 21 ) 


) 


Since- T(t+At) appears on both sides of-the equation, this formulation Is 
termed an imnlicit solution technique-. 

It will be noted that when At approaches Infinity, equation (21) re- 
duces. to the steady state equation discussed in Section 2.3*1 ► Hence, 

If can be deduced (and it has been proven in reference 3) that backward 
differencing is stable for any At. In. addition, backward differencing is 
subject to the-same set- of Iteration termination.crlteria^N and 6, as were- 
applled to the steady state solution. Stability for large At*, however, does 
not preclude the possibility of an oscillating solution*- so a damping 
factor,- c, is also provided to ensure the eventual convergence to the 
solution. 

The great advantage of backward differencing lies in the effective- 
selection of At. During periods of rapidly varying boundary conditions ,. 

At may be reduced to as small a value as desired, and then, during periods 
of slowly^ drifting boundary conditions, It can be enlarged to a compatible 
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vaUie, wt thout regard for-Jthe minimum time constant, t, of the network. 

It alsa happens that when a large At is appropriate,, the nuiriber- of 
Iterations necessary to satisfy the relaxation criterion, 6, Is relatively 
smaU* On the other, hand,, when At IsoompressecLto less than the CSGMIN 
(t), the number of Iterations. per time ste p will still be large when com- 
pared .to forward -differencing* 

2.3.3 Summary of Other Techniques - 

Numerous other approaches to formulating finite, differencing solution* 
are available (reference 5), However, all of these general ly- amount - to 
modifications, or. combination* of basl c forward and backward differencing. ... 
Central differencing, for example, i*.an. implicit technique which computes 
the current-iteration temperatures as the arithmetic average of the forward 
and backward differencing predictions. As another example* exponential 
prediction uses the forward differencing equation but-provldes stability 
for all At by exponentially "derating 1 ' the Q-^ terms according to the 
ratio of to At. 
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2 JL. Modeling Parameters 

The solutions_to the heat transfer equation developed in the pre- 
ceding sections required that continuous variables-.be quantized. Spatial 
variables were quantized through, the artifices of nodes and conductors, 
and time was quantized Into "time steps" at size At. By assuming, for the 
sake, of discussion, that all-nodes are. cubical * with side ax^ tbeo thls 
Ax_cc.n be used as a general, spatial quantum, and_1t_can_be . related to the 
time-quantum,. At^._ For example, the forward-differencing stability 
criterion, t, (and therefore- the. maximum time. step) is related to ax 
as shown in. equation (21): 


one dimensional: m * 2 
two dlmenslonaT: m » 4 
three dimensional: m = 6 


Since the finite differencing solution approaches the exact solution 
as ax and-At approach zero. It Is logical to ask If anything Imposes a 
minimum on these values. The answer is yes: cost and computer core memory 
space. ..Clearly, the latter constralnt-restrlcts ax to -an explicit non- 
zero minimum. -That is, a small Ax means a large, number of nodus and 
conductors, and the computer's memory must" contain enough-space to hold 
all of the parameters (eapacl tance ^temperature^ conductance, etc.J-assa- 
clated with these modeling elements. In addition to using much more-com- 
puter time for analyses, a large model also eosts more to develop than a 
small one. 

The time step can be chosen as small as desired with a consequent 
increase in the computer run time required for analysis. The relationship 
between run time and time step size is linear for forward differencing 
because exactly- one "iteration" Is requl-red for each time step. For 
Implicit methods, however, the relationship is hot as predictable, because 
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the number of-lteratlons Is controlled by- the convergence criteria,. N 
and 6. In addition, the greater the number of iterations, and the greater 
the number of nodes processecLdurlng each Iteration, then the greater 
will be -the susceptibility of the answers to computer round-off error. 

Looking at the problem from- the other direction. Ax may also be chosen 
large enough to Include -the entire thermal system Intone node^ The_t1me 
step, however, at least.for forward differencing. Is limited in size to the 
CSGM11L This .limit Is not Imposed on -the Implicit solution techniques, 
although adamped oscillatory response may result when At is too large. 
Since such, an oscillating response can be critically damped by using the 
damping factor, t. there. is no definitive equation for the maximum At 
that may be used successfully with Implicit solutions. To further compli- 
cate this understanding, there Is no definitive equation for the value 
of e which will yield the most effective damping, although a value of 0.5 
has been routinely used with good results. 

To syitfrtarlze, the thermal math-modeler is faced with the task of 
designing a model and selecting a solution technique which will yield 
good, stable answers for the least cost. To do this, the modeler must 
choose values for the following parameters: 

ax - node size 
At*- time step 
N - Iterations 
6 --convergence criterion 
5 - damping factor 

Forward differencing, which does net require N, 5, or c, and which defines 
at least the maximum At, Is often chosen as the solution technique simply 
because It reduces- the number of parameters which must be juggled about. 

In the hopes of leading to a more logical rationale for selecting model 
and solution parameters, a case stu4y using various methods Is presented 
In Section 3.0. 
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3.0 OPERATIONAL PARAMETER RELATIONSHIPS 


Perhaps the Interrelationship between the lumped parameter system- 
and the- particular solution-techniques that are usually employed- In the 
problem solution can be best displayed with examplesr 

3.1 One-Dimensional Bar of Metal 

The first greup ofc-cases treats. a bar. ef aluminum as shown in- 
FI gyre 2-4 . _ 



Figure 3-1. Bar of Metal 


A boundary temperature of 1000°F Is applled-at x * 0 and_a temperature ef 
0°F Is applied at x ■ 12^. The length of the_bar- radiates to deep space- 
at -46<EF. The bar Is Initially at 0-°F and the transient-response 1s- 
seught at two points: x » 1 and x ® 6. The following cases were run: 


Case No. 

Solution Method- 

# Nodes 

AX 

At 

N 

6 

£ 

1 

Forward 

12 

to 

t 




2 

Forward 

IE 

-1*0 

0.5t 




3 

Forward. 

12 

1.0 - 

Q.It 




4 

Central 

12 

1.0 

T 

50 

to 

1.0 

5 

Centra] - — 

12 

1.0 

3t 

50 

to 

1.0- 

6_~ 

Central 

12 

to 

101--. 

50. 

- - to 

to_ 

7. 

Backward 

12- 

to 


50 

_t0 

ta 

8 

Backward 

12 

1*0 

3t - 

-SO 

— tOL _ 

to 

9 

Backward 

12 

to 

10 T— 

50 

..-.to 

l.O_ . 

10 

Forward- 

3 

4*0 

T 




n 

Forward- 

6 

2.0 

T 




12 

Forward _ 

24 

0.5 

T 




13 

Central 

3 

4*Q 

3f 

50 

1.0 

to 

14 

Central 

6 

2.0 

3t 

50 

to 

to 

15 

Central 

24 

0.5 

3t 

sa 

to 

to 

16 

Backward 

3 

4.0 

3t 

50 

to 

to 

17 

Backward 

6 

2.0 

3t 

50 

1.0- 

to 

18 

Backward 

24 

0.5 

3t 

50 

to 

to 

19 

Backward 

24 

0.5 

T 

50 

to 

to 

20 

Backward 

- 24 

0.5 

t 

50 

0.05 

1*0 


Ignoring gases 19 and 20* theJwlk-of the-xases. ean-he_jd1vided into two 
groups r-{l) time step_yar1atlons» and (Z) node- size variations* Jhe results- 
are shown in Figures 3-2, 3-3*_3-4 v and 3*5.. -These figures plot: accuracy vs 
cost In terms of temperature deviation- vs node-iterations*- Temperature devia- 
tion Is stated-in- degrees F relative to a stated reference value. This value 
represehis-the apparent-correct temperature which was arrived at by corre- 
lating the collection of data points produced by the various cases-. The 
reference temperature is not (necessarily) the correct analytic answer. No 
attempt Was made to find an analytic solution to -the problem. 

The product of the number of hodes times the number of Iterations yields 
a figure that Is closely related to the actual computer cost lnv©lyed-fn 
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analyzing a modef. Fo.* a given computer, the number of node-iterations 
which can be run In a given amount of time is fairly equal regardless of 
the solution technique. Actually, the comparison of node- Iterations for 
a one-d1men$1ortal-model versus node-iterations for a three dimensional 
model is-not strictly reasonable since Jthe latter model will contain more 

conductors per node. However, for^rder*of-magnitude comparisons, node- 

iterations represents a simple and useful measure of cost. 

Figure 3-2 shows -the results for x * 1 inch ancLt »- 25 seconds. The 
effects of -time-step variations and solution techntqye_on a 12 node model 
are showru As. might: be_expected, it appears that the smaller the time 
step, the greater the accuracy- and the greater, the -cost. ...Since, a host-of 
minimax hypotheses could be formulated from this figure, all of which or 
none of which might be true, no further hypotheses will be stated here. 

As Indicated in the Introduction, It remains for the reader to abstract from 
this figure as much meaningful Information as he Is able. 

Figure 3-3 Is similar to Figure 3-2, except that the effects of node 
size variations are emphasized. Since the At and 6 for the implicit routines 
were held constant over the various cases, ftvo additional cases are plotted 
for a smaller Atand fi. Again, as might be expected, it appears that the 
smaller the node size, ax, andJience, the greater the number of nodes, 
then the greater the accuracy and the greater the cost. Further conclu- 
sions will be left for the reader to formulate. 

Figure 3-4 shows the effects of time step variations -and solution 
technique an the results for x s.6 Inches and t * 100-seconds. It will 
br noticed that the temperature dispersion Is net as great as for x * 1 
Inch and-t * 26 seconds.— This is to be expected since the mass of the bar 
from x - 0 to x » 6 serves, as a much more effective low: pass filter, 
which implies greater accuracy fora given sampling frequency (l.e. At}. 

Figure 3-5 is similar to Figure 3-4 except that It highlights the 
effects of node size variations. For a given solution technique, It was 
expected that all points would lie on a smooth curve. This did not hold 


for backward differencing with 24 nodes. . To discover why, two additional 
cases were run with first, the time step,_ar.d then the time step and 
relaxation crlterion-redueed. The results Indicate that -the roll-off at 
large node-iterations was due to the accumulated relaxation error. 
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TEMPERATURE DEVIATION 
(FROM 600 ®F) 



Figure 3-2. ge-Effect of aJfarlatlon In Size of Time Step ffcr Twelve 
Node Model at t«2S seconds and x-ltnch 
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TEMPERATURE DEyiATION 
(FROM 660 °F) 









TEMPERATURE DEVIATION 
(FROM 186 °F) 


Figure 3-4. The Effect of a Venation In Sl2e of Time Step for Twelve 
Node Model at t-ttlG seconds and x-6 Inches 
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NODE- ITERATIONS 







TEMPERATURE DEVIATION 
(FROM 186 °F) 




3.2 Other One- Dimensional Cases 



In the discussion of Section 3J, dimensions were retained so that the - 
reader would- have some Intuitive understanding of the relative magnitudes 
of the various parameters. However, since the basic heat transfer equations 

can be cast In a dimensionless form, the results of- Section 3.1 can.be 

applied to other one-dimensional configurations, as required. 

3.3 . Two-DI mens ional Plate of Metal 

This group of cases treats a plate of aluminum, nine Inches on each 
side and .one inch thicks Boundary conditions are Imposed as 1" me... ted In 
Figure 3-6, and the temperature at the center of the plate, . 100 

seconds. Is desired. 



Figure 3-6. BoundaryjCondltlons on Twe-Dlmenslonal Plate 

The results of a variety of analyses of this plate are presented In 
Figure 3-7 and 3-8. The results In Figure 3-7 appear to be Incorrectly 
plotted since they appear to converge to a point 10 degrees from the 
"assumed" answer. However* the results shown in Figure 3-8 confirm that 
the solution must be elosef to 290°F than 280*F. The difference occurred 
because the model used for the cases displayed In Figure 3-7 contained 
only 36 nodes, ahd this number was evidently Insufficient to accurately 
represent the character of the heat florin the plate. This problem did 
not arise for the one- dimensional cases In Section 3.1 because the heat 

I) 
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X * 

flow, by- definition, always followed the conductor: paths exactly. In 
selecting. finer and f 1 ne r- tnashes f or- the two dimensional- case* it is 
evident -that-the modelled heat flow paths (1 .e^_ the conductors) will lie 
eloser and closer to. the actual paths, with a consequent- increase In 
temperature prediction accuracy. 

3.4 Other Two-Dimensional Cases 

The principles of similarity and non-dimensional analysis may. be. 
applled to the re sultsof. Section 3.3 for other combinations of material 
type and. dl mens 1 ons . 
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TEMPERATURE DEVIATION 
(FROM 290 *F) 
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TEMPERATURE DEVIATION 
(FROM 290 *F) 
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